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INTRODUCTION 


In  practice,  there  often  do  arise  situations  in  which  the 
researcher  wishes  to  choose  the  best  (or  the  t best)  from  a group 
of  k populations,  where  the  "bestness"  of  a population  is  based  on 
ranking  them  according  to  the  value  of  some  characteristic  of 
interest.  In  such  situations,  the  classical  tests  of  homogeneity 
are  inadequate  in  the  sense  that  they  have  not  been  designed  to 
answer  several  possible  questions  in  which  the  researcher  may 
really  be  interested.  Bahadur  [i],  Hosteller  [54]  and  Paulson 
[ 59]  were  among  the  first  research  workers  to  recognize  the 
inadequacy  of  the  usual  tests  of  homogeneity  and  to  consider 
more  meaningful  formulations  in  order  to  answer  these  questions. 
These  developments  set  the  stage  for  the  early  investigations 
of  multiple  decision  problem  which  have  now  come  to  be  known 
as  selection  and  ranking  problems. 

In  the  theory  of  selection  and  ranking  procedures,  there  are 
two  basic  formulations  to  the  problem.  The  first  one  is  called 
the  indifference  zone  formulation  due  to  Bechhofer  [13  ] and  the 
other  is  the  subset  formulation  due  to  Gupta  [29].  The  goal 
of  the  basic  problem  in  the  formulation  of  Bechhofer  is  to  choose 
one  of  the  populations  as  the  best.  The  researcher  is  required 
to  specify  an  "indifference  zone"  in  the  parameter  space  and  the 
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procedure  determines  the  smallest  sample  size  so  that  a certain 
probability  condition  is  satisfied  whenever  the  unknown  parameters 
lie  in  the  "preference  zone".  For  example,  if  we  are  interested 
in  selecting  the  population  with  the  largest  mean  in  N(0.,1), 
i=l,...,k,  the  indifference  zone  is  the  set  of  all  0's  such  that 
the  largest  and  the  second  largest  differ  by  an  amount  £ 6*, 
whereas  the  "preference  zone"  is  the  set  of  all  0's  such  that  the 
above  difference  > S*.  Other  contributions  to  this  indifference 
zone  formulation  have  been  made  by  Bechhofer  and  Sobel  [16], 

Sobel  and  Huyett  [70  ],  Barr  and  Rizvi  [l2],  Desu  and  Sobel  [22], 

Bechhofer,  Kiefer  and  Sobel  [15],  Santner  [64]  and  others.  The 
goal  of  the  basic  problem  in  the  subset  selection  formulation  of 
Gupta  is  to  select  a subset  of  the  given  populations  which  depends 
on  the  outcome  of  the  experiments  and  is  not  fixed  in  advance 
such  that  it  includes  the  best  population  with  a specified  minimum 
probability  regardless  of  the  unknown  configuration  of  parameters, 
i.e.,  over  the  whole  parameter  space.  Some  recent  results  in  the 
area  of  subset  selection  formulation  are  Gnanadesikan  and  Gupta  [28], 

Gupta  and  Studden  [43],  Gupta  and  Panchapakesan  [38], 

Gupta  and  Santner  [41],  Huang  [46]  and  Wong  [72]. 

Many  problems  in  reliability  can  be  considered  in  the  context 
of  selection  and  ranking  problems.  For  example,  one  may  wish  to 
choose  one  or  more  of  the  several  systems  or  components  which  has 
the  largest  mean  life  or  the  largest  median  life.  In  general, 
reliability  problems  also  deal  where  the  distributions  are  unknown 
but  assumed  to  belong  to  a class  of  distributions  such  as  that 
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having  an  increasing  failure  rate  (IFR).  Such  distributions  form 
a special  cases  of  what  are  now  cormonly  known  as  restricted 
families  of  probability  distributions.  The  investigations  of 
Barlow  and  Gupta  [5]  form  the  initial  efforts  on  ranking  problems 
for  such  families.  Other  contributions  to  selection  problem  for 
the  restricted  families  of  probability  distributions  are  Gupta  and 
Panchapakesan  [39,40]  and  Patel  [58].  This  area  of  research  still 
remains  largely  unexplored. 

The  main  investigation  of  this  thesis  is  to  propose  and  study 
selection  procedures  for  some  problems. 

Chapter  I deals  with  some  selection  and  ranking  procedures 
for  restricted  families  of  probability  distributions.  In 
Section  1.1,  definitions  of  various  partial  orderings  on  the 
space  of  distributions  are  given.  In  Section  1.2,  we  propose  and 
study  a selection  rule  for  distributions  which  are  convex-ordered 
with  respect  to  specified  distribution  G.  Some  properties  of  this 
selection  rule  are  discussed.  The  asymptotic  relative  efficiencies 
of  this  rule  with  respect  to  some  other  selection  rules  are  derived. 
Section  1.3  deals  with  the  selecting  the  best  population  using  the 
indifference  zone  approach.  In  Section  1.5,  we  propose  and  study 
a selection  rule  for  distributions  which  are  s-ordered  with 
respect  to  G where  we  are  interested  in  the  scale  parameter  case. 

The  estimation  of  ordered  parameters  from  the  k unknown  distribu- 
tions is  discussed  in  Section  1.8. 

Chapter  II  discusses  some  interval  estimation  problems  from 
k populations.  We  are  interested  in  finding  the  smallest  sample 
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size  N to  be  chosen  from  each  population  such  that  the  probability 

that  a given  confidence  interval  I which  is  based  on  Tr. = max  T. 

l<i<k  ’ 

(T^  is  an  appropriate  statistic  from  i-th  population)  contains  at 
least  one  good  population,  is  at  least  P*.  where  P*  is  a specified 
number,  0 < P*  < 1 . Also  we  are  interested  in  finding  the  smallest 
sample  size  N such  that  the  probability  that  I contains  all  good 
populations  (or  excludes  all  bad  populations)  is  at  least  P*. 

Section  2.2  deals  with  the  above  problems  for  the  location  parameter 
case.  The  infima  of  coverage  probabilities  are  obtained.  For  scale 
parameter  case,  the  above  problems  are  investigated  in  Section  2.3. 
In  Section  2.4,  we  illustrate  the  above  results  by  means  of  two 
examples. 

In  some  situations,  one  has  to  deal  with  a vector-valued 
parameter  associated  X^. . In  such  cases  one  may  consider 
comparing  populations  or  's  in  terms  of  majorization  and  weak 
majorization  of  these  vectors.  Chapter  III  deals  with  such 
problems.  The  parameter  space  is  partially  ordered  by  means  of 
majorization  or  weak  majorization.  Selection  procedures  are 
proposed  and  studied.  In  Section  3.2,  a class  of  procedures 
for  selecting  the  best  population  is  defined.  A sufficient 
condition  is  obtained  for  the  infimum  of  the  probability  of  a 
correct  selection  to  be  Schur-convex  in  A.  Also  another  sufficient 
condition  for  the  same  infimum  of  the  probability  of  a correct 
selection  to  be  nondecreasing  and  Schur-convex  in  A is  obtained 
in  Section  3.3.  Section  3.5  and  3.6  deal  with  selection  procedures 


for  multivariate  normal  distributions  in  terms  of  majorization 
and  weak  majorization.  Various  cases  corresponding  to  the 


known  or  unknown  common  covariance  matrix  s are  studied, 
of  these  selection  procedures  are  also  established. 
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CHAPTER  I 

SELECTION  PROCEDURES  FOR  RESTRICTED  FAMILIES 
OF  PROBABILITY  DISTRIBUTIONS 

1.1  Introduction 

In  many  problems,  especially  those  in  reliability  theory,  one 
is  interested  in  using  a model  for  life  length  distribution  which 
belongs,  for  example,  to  a family  of  distributions  having  increasing 
failure  rate  (IFR),  or  increasing  failure  rate  on  the  average  (IFRA), 
Such  distributions  form  special  cases  of  what  are  now  commonly  known 
as  restricted  families  of  probability  distributions.  These  are 
defined  more  precisely  later  in  this  section.  The  idea  of  using 
such  families  stems  from  the  fact  that  in  many  cases  the  experimenter 
cannot  specify  the  model  (distribution)  exactly  but  is  able  to  say 
whether  it  comes  from  a family  of  distributions  such  as  IFR,  IFRA. 
Families  of  probability  distributions  of  these  types  have  been 
studied  by  several  authors,  see,  for  example,  Barlow,  Marshall  and 
Proschan  [7],  These  authors  have  mainly  concerned  themselves  with 
probabilistic  aspects  of  these  distributions.  To  some  extent,  there 
have  been  some  investigations  dealing  with  statistical  inference  for 
some  of  these  families;  see  for  example,  Barlow  and  Proschan  [8],  [9], 
Barlow  and  Doksum  [3]. 
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In  this  chapter  we  are  interested  in  studying  multiple  decision 
procedures  for  k (k  2)  populations  which  are  themselves  unknown  but 
which  are  assumed  to  belong  to  a restricted  family.  First  we  give 
some  notations  and  definitions.  A binary  ordering  relation  (-<)  is 
called  a partial  ordering  in  the  space  of  probability  distributions 
if 

(a)  F -<  F for  all  distributions  F,  and 

(b)  F <G,  G imply  F <H. 

Note  that  F -<  G and  G < F do  not  necessarily  imply  F ^ G. 

We  now  define  some  of  the  special  order  relations  of  interest  to 
us  (see  Barlow  and  Gupta  [5]). 

(i)  F is  said  to  be  convex  with  respect  to  G (written  F < G)  if 

c 

and  only  if  G'V(x)  is  convex  on  the  support  of  F. 

(ii)  F is  said  to  be  star-shaped  with  respect  to  G (written  F < G) 

★ 

if  and  only  if  F(0)  = G(0)  = 0 and  ^ is  increasing  in 
X ^ 0 on  the  support  of  F. 

(iii)  F is  said  to  be  r-ordered  with  respect  to  G (F  < G)  if  and 

1 G”^F(x1  ^ 

only  if  F(0)  = G(0)  = ^ and  ^ is  increasing  (decreasing) 
for  X positive  (negative)  on  the  support  of  F. 

(iv)  F is  said  to  be  s-ordered  with  respect  to  G (F  < G)  if  and 
only  if  F(0)  = G(0)  = ^ and  G’^F  is  concave-convex  about  the 
origin,  on  the  support  of  F. 

If  G(x)  = 1-e"^^,  X ^ 0,  then  F < G is  equivalent  to  saying  that 
F has  increasing  failure  rate  (IFR).  The  class  of  IFR  distributions 
has  been  studied  by  Barlow,  Marshall  and  Proschan  [7]. 
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Again  if  G(x)  = 1-e  , x ^ 0,  F •<  G is  equivalent  to  saying 

★ 

that  F has  increasing  failure  rate  on  average  (IFRA).  The  class  of 
IFRA  distributions  has  been  studied  by  Barlow,  Esary  and  Marshall  [4]. 
The  r-ordering  has  been  defined  and  investigated  by  Lawrence  [48], 

The  s-ordering  and  c-ordering  have  been  studied  by  Van  Zwet  [71  ] 
and  Lawrence  [48]. 

In  the  statistical  literature,  selection  problems  for  restricted 
families  were  first  investigated  by  Barlow  and  Gupta  [5].  Some 
further  results  in  this  direction  and  a review  of  some  important 
results  concerning  inequalities  for  restricted  families  and  problems 
of  inference  for  such  families  have  been  given  by  Gupta  and 
Panchapakesan  [39].  The  selection  of  the  population  with  largest 
a-quantile  from  distributions  which  are  star-shaped  with  respect  to 
the  folded  normal  distribution  has  been  considered  by  Gupta  and 
Panchapakesan  [40].  In  a recent  paper,  Patel  [58]  has  studied  the 
selection  of  IFR  populations  which  differ  only  in  the  scale 
parameters.  His  procedure  is  based  on  the  total  life  statistic 
until  r-th  failure. 

In  Section  1.2,  we  propose  and  study  a selection  rule  for 
distributions  which  are  ordered  with  respect  to  a specified 
distribution  G assuming  there  exists  a best  one.  Some  properties 
of  this  rule  are  discussed.  The  infimum  of  the  probability  of 
a correct  selection  is  obtained  and  an  asymptotic  expression  is 
also  given.  We  also  study  the  asymptotic  relative  efficiencies  of 
this  rule  with  respect  to  some  selection  procedures.  Section  1.3 


deals  with  the  selecting  the  best  population  in  the  frame  work  of 
Section  1.2  using  the  indifference  zone  approach.  In  Section  1.4, 
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we  propose  a selection  procedure  for  distributions  that  are 
ordered  with  respect  to  G.  Section  1.5  deals  with  selection  proce- 
dure with  respect  to  the  means  for  distributions  that  are  s-ordered 
with  respect  to  G where  we  are  interested  in  the  scale  parameter 
case.  The  distribution  of  V.  (see  Section  1.5)  is  also  investi- 
gated. Section  1.8  deals  with  estimation  of  ordered  parameters 
from  k unknown  distributions  where  we  are  interested  in  the  scale 
parameter  case. 

1.2  Selection  rules  for  distributions  ordered  with  respect  to 
a specified  distribution  6. 

Let  5 be  the  class  of  absolutely  continuous  distribution 
functions  F on  R with  positive  and  right-{or  left-)  continuous 
density  f on  the  interval  where  0 < F ' 1.  It  follows  that  the 
inverse  function  F’^  is  uniquely  determined  on  (0,1).  We  take 
F’^(O)  and  F'^(l)  to  be  equal  to  the  left  hand  and  right  hand 
endpoints  of  the  support  of  F.  For  F,  G € ir,  consider  the  following 
transformation  (see  Barlow  and  Doksum  [ 3]) 

1 1 

(1.2.1)  H:'(t)  = / g[G‘T(u)]du,  0<t<l. 

F’^0) 

We  assume  that  G is  always  fixed.  Since  Hp^  (the  inverse  of  H)  is 
strictly  increasing  on  [0,1],  Hp  is  a distribution.  We  know  that 
F < G if  and  only  if  Hp  is  convex  on  the  interval  where  0 < Hp  < 1. 
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Since  G is  assumed  known  we  can  estimate  Hp  by  substituting  the 
empirical  distribution  of  F;  that  is 

F;'(t) 

(1.2.2)  H''(t)  = = / 9[G‘'F^(u)]du 

" F’\o) 


(1.2.3)  H-'(t)  = /"’"  g[G-'F„(u)]du  = .1  g[G-'(i^)](X,_„-X,.,^„) 

F''(0) 

where  X.  is  the  i-th  order  statistic  in  a sample  of  size  n from 
1 » n 

F and  X^,  „ 0. 

0,n 

If  G(x)  = l-e’’^  for  x ^ 0,  then  (1.2.3)  can  be  written  as 


say  that  X^  ^ +...+  X^_-j  ^ + (n-r+l)X^  ^ is  the  total  life  statis- 


tic until  r-th  failure  from  F.  Thus,  H (^)  = -I-  { 


-l/Tx  _ 1 .total  life  statistic, 


n 'n'  n until  r-th  failure 


(A)  Selection  procedure  and  its  properties 


Let  TT.| , . . . ,7T|^  be  k populations.  The  random  variable  X.j  associa- 
ted with  iT^.  has  distribution  function  F.j , i = l,2,...,k,  where  F^.  e 3 
(i  = l,...,k).  Let  Fj-|^j  denote  the  cumulative  distribution  function 
(c.d.f.)  of  the  "best"  population.  For  example,  if  we  are  interested 
in  the  quantile  selection  problem,  the  "best"  population  may  be 
defined  as  that  Fj-j^j  (unknown)  which  satisfies  (a)  below. 

We  assume  that 

(a)  ^ l,...,k-l; 


n 


(b)  there  exists  a distribution  G such  that 


where  < denotes  a partial  ordering  relation  in  the  space  of 
probability  distributions.  It  should  be  pointed  out  that  the 
condition  (a)  above  may  also  imply  that  is  the  distribution 

with  the  largest  (smallest)  parameter.  For  example,  if 

^i(^)  ^ 2 0,  0^  > 0 (i  = l,...,k),  then  F|-|^j  is  the 

distribution  with  the  largest  o^-.  We  are  given  a sample  of  size  n 
from  each  (i  = l,...,k).  Our  goal  is  to  select  a subset  from 
the  k populations  so  as  to  include  the  population  with  Let 

. =^F=(F^  . ,F|^) ; s aj  such  that  F^(x)  ^ 


(1.2.5) 

(1.2.6) 


T.  = 


1 for  i = l,...,k. 


T = 2 a.  Y. 


j = l 


J J ,n 


where  X.  . is  the  i-th  order  statistic  from  F.,  Y.  is  the  j-th 
i;j,n  i’  j,n 

order  statistic  from  G,  r is  a fixed  positive  integer  (1  r ^ n), 
= g G'^(^)  - g G'^(i)  for  j = l,...,r-l 
and  = g G'^(^). 

For  selecting  a subset  containing  we  propose  the  selection 

rule  R.|  as  follows: 

R-| : Select  population  if  and  only  if 
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(1.2.7)  ^ ■ 

where  c.|  = c^(k,P*,n,r)  is  some  number  between  0 and  1 which  is 
determined  as  to  satisfy  the  probability  requirement 

(1.2.8)  inf  P1CS;R^  ! > 

where  CS  stands  for  a correct  selection,  i.e.,  the  selection  of 
any  subset  which  contains  the  population  with  distribution 
Let  be  associated  with  and  W^(x)  be  the  c.d.f.  of 
from  We  now  state  a few  preliminary  lemmas. 

Lemma  1.2.1.  (Lehman  [49]  p.  112)  Let  F(x)  be  a distribution 
function  on  the  real  line.  If  v(x)  is  any  nondecreasing  function 
of  X,  then  E ij'(x)  is  a nonincreasing  function  of  F,  i.e.,  if  F^(x) 
F^(x)  for  all  X,  then  /ii'(x)dp2(x)  < /i;;(x)dF^  (x). 

Lemma  1.2.2.  Let  Xp...,X^  be  i.i.d.  with  distribution  F(x).  Let 
; be  a function  of  Xp...,X^  which  is  nondecreasing  in  each  of  its 
arguments.  Then  Ei^(X^ , . . . ,X^)  is  nonincreasing  function  of  F. 


Lemma  1.2.3.  (Gupta  and  McDonald  [35])  Let  X = (X,n,...,X,  ,..., 

- ii|^  in^ 

X.  X.  ) be  a vector  valued  random  variable  of  )'  n.  (-  1) 

kl  kn,^  1 - 

independent  components  with  X^.j  having  the  distribution  F^(x), 


j = l,...,n^,  i = l,...,k.  Let  i|i  be  a function  of  X.|p...,X.|^  ,..., 

Xkr'’‘,X|<n  fixed  i,  is  a nondecreasing  (nonincreas- 

k 

ing)  function  of  X^. .j , . . . ,X . ^ when  the  other  components  of  X are  held 
fixed.  Then  Eii(X)  is  a nonincreasing  (nondecreasing)  function  of 


F.  . 
1 
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The  following  two  lemnias  are  from  Barlow  and  Proschan  [8  ]. 
n n 

Lemma  1.2.4.  (a)  J a,Xj  - 1(0)  • J a-  [-{-(Xj  - ^*(0)]  for  all 

■"  ^ = 1 1 1 - 1 1 

0 • X-|  ^ b and  for  all  convex  4.  on  [0,b]  if  and  only  if 

n 

0 ^ A.  ■ I a . V 1 for  i = l,2,...,n. 
j = i ^ 

(b)  ,’>()'  a. X. )-  4(0)  ; a.[.j(X. ) - ; (0)]  for  all  0 £ X,  -. . .1  X 

i=l  1 1 i=l  1 1 _ _ ' ' 

and  for  all  convex  4.  on  if  and  only  if  A-j  ^1,  A2  1 1 , . • . , 

^ ’ "^k+l  - 0,. . . .>A^  ^ 0 for  some  k (0  k ;<  n) . 


Lemma  1.2.5.  (a)  4'(  J a-X.)  } a.  4,(X.) 

i=l  ’ ’ “ i=l  ’ ’ 

for  all  star-shaped  4>  on  [0,b]  and  all  0 ;;  X.j  X^  ^ b if  and 

only  if  there  exists  k (1  ^ k ^ n)  such  that 

0 T:  1-  • -1  £ 1 and  Aj^^^  =. . .=  A^  = 0. 

n n 

(b)  4(y  a-X.)>  y a.4>(X.)  for  all  0 £ X-,  £•  • -f.  X and  for  all 
i=l  ^ “ i=l  ’ ^ I - n 

star-shaped  4 on  if  and  only  if  there  exists  k (1  £ k £ n) 


such  that 


£•••!  ^k+1  "•••"  ^n  " °- 


Lemma  1.2.6.  Let  F-| , be  two  distribution  functions  such  that 
F^(x)  > F2(x)  Vx. 

Let  T.  = ):  b.  X . i = 1,2 

1 J i>J,n 

where  b . > 0 for  j L A,  a <_  { 1 ,2, . . . ,n } and  X.  . is  the  j-th 
J ‘ » J > *■ 

order  statistic  from  F^,  i = 1,2. 


P[T^  £ x]  £ P[T2  : x]. 


Proof. 


f 
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Let  (X.^.-.-.X,^)  = 


0 


if  T.  ^ X 


otherwi se 


where  X.i,...,X^^  are  n observations  from  (i  = 1,2). 

Since  (X^- ^ • .X^^)  is  nondecreasing  in  each  of  its 
arguments,  by  Lemma  1.2.2  we  have 

F ’vl)i(X^1,...,X^^)  ^ L i(<2(X2i»---,X2^) 

That  is  P[T-|  i x]  ^ P[l2  ^ x] 

This  proves  the  Lemma. 

We  now  state  and  prove  following  theorem  which  is  more  general  than 
that  of  Patel  [58]. 

Theorem  1.2.1.  LetF^,  G € ir,  F^(x)  i Fj-j^j(x)  Vx  and  i = l,2,...,k, 
Fr, -i(O)  = 0 and  Fr, . < G. 


[k] 


[k] 


If  3j  0 for  j = 1,2,. ..,r,  G ^(0)  £ 0,  g G"^(0)  £ 1 and  a^  ^ c-j 
then 
(1.2.9) 


P[CS|R^]  >/  G!f’’'(|-)dG.^(x) 

G’^(O)  ^ 


where  Gy(x)  is  the  c.d.f.  of  T, 


Proof.  P[CS|R^]  = P[T(^)  > c^  l^.y  i ^ k] 

k-1  k 1 X 

= / II  )dw.(x) 

0 i = l ^ ‘^l 

> / w!^{f-)dW.  (x)  (By  Lemma  1.2.6) 
■ 0 1 

= P[z,  > c^Zj,  j ^ k] 

where  Zi,...,Z|^  are  i.i.d.  with  c.d.f.  W,,(x). 


1 
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Let  cf(x)  = G' 


Ncte  that 
(1.2.10) 


^ y a.  Xt  . 

1 St  jii  J i;j,n 


i ~ l(...iki 


where  • is  the  j-th  order  statistic  in  a sample  of  size  n 
from  i = 1 , . . . ,k. 

(1.2.11)  P[2.  c,  max  Z.]  = P[cp(^  Z.  ) ^ cf(Z.),  i = l,...,k-l 

llJlk  ^ 1 

r 1 

Since  [ a . = g G (0)  ^ 1 and  a • ^ 0 V j = 1 , . . . ,r , then,  by 
j=l 

Lemma  1.2.4  (a)  and  (1,2.10), 

0.2,12) 

1 1 *" 

Since  — a > 1 and  — J a.  > 1 for  i = l,.,.,r,  we  have,  by 
^1  '^l  j=i  ^ “ 

Lemma  1.2.4  (b)  and  (1.2,10), 

0.2.13)  <^XJ.j_„). 

where  Y.  . is  the  j-th  order  statistic  from  G,  i = l,2,...,k. 

1 » J *n 

Thus  from  (1.2.11),  (1.2.12),  (1.2.13)  and  (1.2.14), 

P[Z^  > c,  Z,]  > P[ a.V|^,_j_„  > c,  ajY,,j 


^ LI 

/ G:‘Uf-)dG,(x) 

G-l(O)  ^ 


This  completes  the  proof. 
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- The  constant  = c-j (k,P*,n,r)  satisfying  (1.2.8)  is  determined  by 

I / a!f'^{^)dG^(x)  > P*  and 

G‘^(0)  ^ S 

g i j 

We  now  consider  two  specific  distributions  G(x).  If  G(x)  = 1-e"^, 

X ^ 0-,  then  we  have  following  result  which  generalizes  the  result  of 
Patel  [58]. 

Corollary  1.2.1.  If  (x)  ^ F^l^j(x)  V x and  i = 1 , . . . ,k,  F^|^j(0)=0, 

F|-|^j  -c  G,  G(x)  = 1-e"^,  X ^ 0 and  n > max{r,  then 

00 

(1.2.15)  inf  P[CS1R,]  = / H'^‘^(MdH(x) 

0 ‘^1 

2 

where  H(x)  is  the  c.d.f.  of  a x random  variable  with  2r  d.f. 

Proof.  If  G(x)  = 1-e'^  then  a.  = i for  j = l,2,...,r-l 

vl  * * 

and  a^  = 1 (n-r+1). 

I 

Also  — a > 1 iff  n > ^ — . 

Cl  r - - 

2 

By  Theorem  1.2.1  and  the  fact  that  2nT  is  distributed  as  x with  2r 

d.f.,  the  result  follows.  ^ 

If  G(x)  = x for  0 < X < 1,  then  we  have  the  following  result  which  is 

a special  case  of  Theorem  2.1  of  Barlow  and  Gupta  [5]. 

Corollary  1 .2.2.  If  F^(x)  ^ F|-|^j(x)  V x and  i = l,...,k,  F|-|^^(0)=0, 

-<  G and  G(x)  = x for  0 < x < 1, 


then 
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(1.2.16)  inf  P[CS|R,]  = fx) 

' 0 1 

where  H^(x)  is  the  c.d.f.  of  the  r-th  order  statistic  from  G. 

Proof.  If  G(x)  = X 0 £ X < 1 then  = 0 for  j = 1 , , . . ,r-l  and 
■ 1 • By  Theorem  1.2.1,  the  result  follows. 

We  state  and  prove  the  following  theorem  about  the  asymptotic  evalu- 
ations of  the  probability  of  a correct  selection. 


Theorem  1.2.2.  If  G c for  all  i = l,..,k  and 

(i)  F.(x)  > Fj-|^j(x)  V X,  i = l,...,k,  and  Fj-j^^  < G, 

(ii)  G(x)  has  a differentiable  density  g in  a neighborhood  of 
its  a-quantile  and  g('i^)  f 0, 

(iii)  g G"^  is  uniformly  continuous  on  [0,1),  G'^{x)  is  convex  and 
there  exists  an  ri,  0 ■ ti  1,  such  that  for  m £ y < 1 > g G ^ (y)  is 
nonincreasing  and  is  nondecreasing  in  y,  then  as  n ^ " 

(1.2.17)  P[CS|R,]  > / g(n  )(— )^]d«>(x) 

I 1. 1 L 1 U 01  ~ 

-or>  I I am 

where  ^ a as  n li  = l-<  and  i>(x)  is  the  standard  normal  c.d.f. 
Proof.  See  the  proof  of  Theorem  1.2.1,  we  have 


(1.2.18)  P[CSiRj  > P[Z.  • c,  max  Z,] 

' 'lrj<k  J 

where  Z.|,...,Z|^  are  i.i.d.  with  c.d.f.  Wj^(x)  and  W|^(x)  is  the  c.d.f. 
of  T 


,2 


(k)- 

By  Theorem  2.2  of  Barlow  and  Van  7wet  [11]  and  (iii). 


18 


(1.2.19)  sup]/  g[G"  Fj-|^^^(u)]du  - / g[G’  F|-|^j(u)]du|-^  a.s. 

x>0  r"l  / n\  r"l  I n\ 


where  is  the  empirical  distribution  of  Fj-j^j. 

Then  we  have  (see  Barlow  and  Doksum  [3]),  for  n large. 


(1.2.20) 


^i  St  ^i ;r,n 


where  Y.  is  the  r-th  order  statistic  from  Hr  and  Hr  (the 

"^[k]  '^[k] 

inverse  of  Hr  ) is  defined  in  (1.2.1).  Since  Fr.-,  < G,  therefore 
F[k]  LkJ  c 

Hr  is  convex.  Since  G’^(x)  is  increasing  and  convex,  therefore 
'^[k] 

G'^Hr  (x)  is  convex.  Since  Hr  •<  G and  G*^(0)  < 0,  therefore 
■^[k]  "^[k]  c 

Hr  ■<  G.  In  a manner  similar  to  the  Theorem  2.1  of  Barlow  and 
■^[k]  * 

Gupta  [5],  we  have 

"[’'k;r,n  i '’t''k;r,„  i ^ 

where  Yf  is  the  r-th  order  statistic  from  G,  i = l,...,k.  From 
1 ;r  ,n 

(1.2.18),  (1.2.20),  (1.2.21)  and  using  the  fact  that 


the  theorem  follows. 

Consistent  with  the  basic  probability  requirement,  we  would  like 
the  size  of  the  selected  subset  to  be  small.  Let  S be  the  size 
of  the  selected  subset.  One  criterion  of  the  efficiency  of 
the  procedure  is  the  expected  value  of  the  size  of  the  subset. 
If  in  addition  to  the  assumptions  of  the  Theorem  1.2.1, 
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we  assume  that  (a)  FrT,(x)  f^rn(r)  all  x > 0, 


i = 1 ,2,. .. ,k,  where 


[1]^M  J-  i 

1 is  given  and  (b)  G < F 


[!]• 


We  have  the  following  theorem. 


Theorem  1.2.3. 

fX> 

(1.2.22)  E(SlRJ  < k / x)dG,(x) 

G '(0)  ' 

where  Gy(x)  is  the  c.d.f.  of  T. 

k 

Proof.  E(S|R,)  = I P[T.  > c,T.,  for  j i]. 
i = l ^ J 

By  Lemma  1.2.6,  we  have 

P[T^  1 Mil  P[T*  > c^TJ,  j = l,..,k-l] 

^k 

where  Tj  has  c.d.f.  W.|(x)  for  j = l,...,k-l  and  ~ has  c.d.f.  W.j(x). 

So  that  E(S|R^)  < k P[T*  > c^Tj,  j = l,...,k-l] 

= k P[Z^5  > c^Zj,  j = l,...,k-l] 

where  Z^,...,Zj^  are  i.i.d.  with  cdf  W.|(x). 

Hence 

E(SiR^)  < k P[^^Z^  > Zj,  j = 1,2,. ...k] 

Let  q(x)  = G'Vpj(x) 

since  G <Fr,T  then  - cf(x)  is  convex, 
c ^ 

Using  the  same  approach  as  in  Theorem  1.2.1 
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and  ~ a > ^ a^,  we  have 
C]  r - r 

Z,^  > Zj.  j = l,2,....k]  ' P[^y  Z*  > Z^,  j = l,....k] 

CO 

-!  Glf-'(^x)dG^(x) 

G‘^0)  ^ 


where  Z|,...,Zj^  are  i.i.d.  random  variables  with  c.d.f.  Gy(x). 
This  proves  the  theorem. 


Before  discussing  properties  of  the  selection  rule  R.j , we  give  some 
preliminary  definitions.  Let  u denote  the  set  of  all  k-tuples 
F = (Fp...,F|^).  Let  Pp(i)  denote 


Pp('i)  = Pp^^(i)  selected [R] 

when  is  associated  with 

Defi ni tion  1.2.1.  A rule  R is  strongly  monotone  in  if 

l' 

Pp(i)  is  ^[i]  other  components  of  F 

\ are  fixed 

( t in  other  components  of 

F are  fixed. 

Definition  1.2.2.  A rule  R is  monotone  means 


Pp(i)  j;!  Pp(j)  for  all  F e with  Fj-^j(x)  ^ • 
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Definition  1.2.3.  A rule  R is  unbiased  if  Pp(i)  £ Pp(k)  for  all  F € n 
with  P[-ij(x)  f[;|(j(x). 

Definition  1.2.4.  A rule  R is  consistent  with  respect  to  ii'  means 
inf  P[CSjR]  1 as  n * «. 

Ti' 

Theorem  1.2.4.  If  a^.  ^0  for  i = l,...,r,  then  R-j  is  strongly 
monotone  in 

Proof. 


Let  i>(x)  = 


1 if  T 


(i)  ^ '(j) 


0 otherwise 


since  a-  ^ 0,  is  nondecreasing  in  each  of  its  arguments  for 
i = l,...,k. 

So  that  i^(x)  is  nondecreasing  function  of  )i » • • • »^(i  )n 
other  components  of  X are  held  fixed  and  is  nonincreasing 
function  of  j )i  ’ • • • j )n  ^ ’ ^ other  components  of  X 

are  held  fixed  where  j )]’••■ ’^(j )n  observations  from 

F[j](j  = l,..,k).  Since  E4>(X)  = Pp(i),  by  Lemma  1.2.3,  it  follows 
that  R.|  is  strongly  monotone  in 

Remark  1.2.1.  (1)  If  a rule  R is  strongly  monotone  in  j for 

all  i = l,...,k,  then  R is  montone  and 

Inf  P[CS|R]  = Inf  P[CSlR] 

“0 

where  = (F  = (F^ , . . , ,F|^) ; F^  =...=  Fj^). 

(2)  If  R is  monotone,  then  it  is  unbiased. 

(3)  If  F^(x)  = F(x,0.j),  i = l,...,k  and  T.  is  a consistent  estimator  of 
e.,  then  R^  is  consistent. 

(4)  If  F^ , G L j,  F.j  ^ G,  i = 1,...,k  and  the  condition  (iii)  (excluding 
G~^(x)  is  convex)  of  Theorem  1.2.2  is  satisfied,  we  can  show  that 
R^  is  consistent. 
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The  selection  of  the  population  with  largest  (i  = l,..,k)  can  be 
handled  analogously.  We  assume  F[-^j(x)  iFpj(x),  i = l,...,k,  and 

Fj-^j  < G.  The  rule  for  selecting  the  population  with  Fj-^j  is 

R2:  Select  population  if  and  only  if 


(1.2.23) 


Co  T-  •'  mm  T. 
^ 1 ~ 1 • I ] 

I'j'k 


where  C2(0  C2:_  1)  is  determined  so  as  to  satisfy  the  basic 

requirement.  In  a manner  similar  to  the  proof  of  Theorem  1.2.1, 
if  F^. , G € ii,  Fj-.j(x)  1 Fj-^j(x)  V X and  i = l,...,k,  Fpj(O)  = 0 
and  Fp  j -c  G and  if  ^0  for  j = 1 , . . . ,r,  G~^  (0)  ^ 0 , g G"^  (0)  ^1 


and  a^  > C2>  then 


(1.2.24)  P[CS!R2]i/  g'^‘^C2x)dG^(x) 

G‘^(0) 


where  G^(x)  = 1-G.j.(x). 

( B ) Efficiency  of  procedure  R -j  under  slippage  configuration. 

We  consider  slippage  configuration  Fj-^.j(x)  = F(j),  i = l,2,..,k-l, 
and  i^[-|<.](x)  = F(x),  0 ■ i ■-  1.  Let  E(S|R)  denote  the  expected  subset 
size  using  the  rule  R.  Then  E(SjR)-P[CS|R]  is  the  expected  number 
of  non-best  populations  included  in  the  selected  subset.  For 
a given  C .•  0,  let  nj^(£)  be  the  asymptotic  sample  size  for  which 
E( S I R) -PICS I R}  = C.  We  define  the  asymptotic  relative  efficiency 
A R E(R,R*,-,)  of  R relative  to  R*  to  be  the  limit  as  c ' 0 of  the 


IL 
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nn(t)  n„(c) 

ratio 7— Y are  (R,R*;^)  = i im  r — r?r- 

nR.TJ 

Under  the  slippage  configuration,  for  large  n,  we  have  (see 
Theorem  1.2.3) 


(1.2.25)  E(S|R^)  = P[CS1R^]  + (k-l)P[T^^j  max  T^.j] 

(1.2.26)  P[Tn\  ' c max  T,.  J - P[Y,  - c,  max  Y-] 

U1  - I --  I 1 

where  Y^,...,Y|^  are  independent  and  Y^.  is  the  r-th  order  statistic 

from  Hr  for  i = l,...,k.  The  right-hand  side  of  (1.2.26)  is 
'^['1 

equal  to 

(1.2.27)  / T(—  - a^h(aj(l  - ^)(-^)^). 


- a,  h(a,  )(1  - ^)(-^)^)d4(x) 
1 I aa 


where  c^  is  the  constant  used  in  defining  R^ , a^  is  the  (unique) 

.«-quantile  of  Hf  (x)  and  h(x)  is  the  density  function  of  Hp  (x). 

^[k]  ^[k] 

From  now  on,  it  is  assumed  that  k = 2.  Therefore  (1.2.25)  reduces  to 


(1.2.28)  E(S|R,)-P[CS|RJ  = T(-h(a  )a  (1  - ^)(-^F(l  + ^)“" 

I I Ot  Ot  ^ 1 ^ 

I aa  C.| 

Let  +0-cJn  q(n  ) k-¥)dt(x)  = P*. 

-TO  I I aa 

Barlow  and  Gupta  [5]  have  proved  that 

, _ 1 2«D  ^ D 3 + n/l  i 

- X * n - :V2  jri  ‘><7> 


where  D 


- T~''(P*)(aa)^ 
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Now,  setting  the  right  hand  side  of  (1.2.28)  equal  to  C and  using 


2^ 


c-|  , we  obtain 


(1.2.29)  n„  (6)  ..  [-(...i)=v-l(C)(l+c^)r+/2  D6a  h(a  )]^- 


[a.^h^(a 


Comparison  with  Barlow-Gupta  Procedure 

Barlow  and  Gupta  [5]  propose  a procedure  R^,  for  a quantile 
selection  problem. 

R2;  Select  population  if  and  only  if 


(1.2.30) 


T . > C3  max  T 
l"j<k 


where  c.  (0  < c ^ 1)  is  chosen  to  sat'sfy  P[CS|R]  > P*  and  T .is 
j — r , 1 

the  r-th  order  statistic  from  where  r £ (n+l)a  < r+1 . A similar 


expression  for  Op  ( e)  (see  Barlow  and  Gupta  [5])  is 
^3 


2^^ 


Hr  (€)  ^ [-(au)'^«''(e)  (1+.;^)‘^  + v^  D-V  f(r  )]^[c^f^(£  )(l-6)^] 

f'  ^ 'Jt  CX  Cl  Cl 

where  f is  the  density  of  F with  unique  a-quantile,  c • 


,2.-1 


(1.2.31)  ARE(R,,R2;^)  = lim  = - 4- 

' a^  h^(a  ) 


If  G(x)  = 1-e  , X - 0 and  F(x)  = 1-e  , x > 0,  then 

2 2 

(1.2.32)  ARE(R,,R2;-)  = < 1 . 

t 

_ 1 
* 2- 


= 0,4803, 


•1 
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Comp a r1 son  with  Gupta  Pro cedur e 

Gupta  [30]  gave  a selection  procedure  for  gamma  populations 

X 

■I  . 's  with  densities  — ^ e ’ x • 0,  o . 0,  i = 1 ,2, . . . ,k. 

r(a)o*:  ’ 

The  procedure  is 

R^:  Select  population  if  and  only  if 

(1.2.33)  X.  > c,  max  X- 

where  is  the  sample  mean  of  size  n from  and  c^  is  the  largest 
constant  (0  <c^  1)  chosen  so  that  P[CS|R^]  ^ P*. 

For  k = 2,  ipj  = • and  = 1 , Barlow  and  Gupta  [5]  have  prove  that 

Or  ( f.) 

(1.2.34)  ARE(R.,R,;  ) = lim 

^ C .0 

a(loq  ■ )^  ii(1  + 

2(1-  )^['-  f(^  )]^ 

■t  X 

It  is  easy  to  show  that 

(1.2.35)  ARE(R^,R^;0  = ARE(R^  ,R3;.- )ARE(R3,R^;f ) 

; log  ^ /x~t  /l+7 j 2 
\'Z  (l--)a  h(a  ) 

t I 

If  G(x)  = 1-e'^  for  x - 0 and  a = 1, 
then 

(1.2.36)  ARE(R,,R.;.S)  = . 

‘ ^ 2(l-s)  a 


Therefore  as  s ► 1 , we  get 


(1.2.37)  ARE(R^,R^;6tl)  = 


It  is  easy  to  see  that 


f'  ’’  2 


(1.2.38)  ARE'R,,R.;^il)i=  1, 


‘ 2 


h.  1,  .<1 


Comparison  of  R ^ and  R ^ from  uniform  distribution 

Suppose  and  two  independent  uniform  populations  with 

distribution  functions  (i  = 1,2). 


X < 0 


F,(x)  = ^ 


0 < X ^ 0. 


X > u . 

1 


where  <•  Oj-^j  = 1. 

A sample  size  n is  drawn  from  each  of  the  two  populations.  Let  T.* 
be  the  total  life  statistic  until  r-th  failure  from  7.  (i  = 1,2) 
where  r (n+1  )-i  • r+1 . 

The  procedure  R^  is  given  by 
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R^:  Select  population  if  and  only  if 


(1.2.39) 


Tt  > Cn  max  T* 
' ' ^ ]<j<k  J 


where  is  chosen  so  that  PCCSIR^]  > P*.  Let  be  associated 
with 

(1.2.40)  E(S1R5)  - PCCSlR^]  = P[Tf^)  > c^TJ^)]  = P[Tj  > / T'] 

where  T^,  are  two  independent  total  life  statistic  until  r-th 
failure  from  uniform  distribution  over  (0,1).  By  Gupta  and 
Sobel  [42], 


(1.2.41) 


N(0,1) 


as  n 


where  u = ,,  u’  = 

2 an-a)n^(2-a)^n^ 

' ^ 4n3  12n2 

= A n where  A = ^ . 

Hence  - « — = B/n  where  B = . 

0 0 2/K 

From  (1 .2.40) , we  have 


T]-u'  Ct-  U-u’  c 


ECSIRj)-  PCCSIRj]  • i ^ (- 

where  are  i.i.d.  with  N(0,1). 


lo-U  Lr  I 

,r  (V-)*<r 


Hence 


E(S|R5)  - PCCSIR^]  = / x-(l-  ^ )B/n'  ]d4>(x) 

5 5 


1 


- (1  - 


/h{M2 

Let  E(SjRg)  - P[CS[Rg]  = t > 0,  we  obtain 


n,2,«)  (L.i,  4; 

5 / cg 


Note  that 


inf  PCCSiRg]  = P[T|  > Tp 
where  T-j,  are  defined  as  above. 

P[Tj>cT’].p[!l^,C5t^l*(VUfl 

PCZ-,  i CgZ^  + (Cg-l  )B/fn 

where  Z-j , I2  are  i.i.d.  with  N(0,1). 

Hence 


P[Tj  > cT']  = / t(^  X -(1  - ]-)B/n)d4>(x) 


-0  - ™)B,^ 

Cr 

= •[— — i ] 

5 

Setting  inf  P[C5|R^]  = P*,  we  obtain 

- (1  - (P*)/mi 

^ ^5 


■ Vp*1 


(l-c^)r^  = d/i  + cl  Where  0 = 


We  see  that  p-  1- 


rr  1 + 


/2  D 
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From  (1.2.42), 


(1 


/FT 


/ii 


-1, 


/Ff(l  - i)  + {^+  1}^ 

(S 


Thus 


0.2..:<'  n^U)  . »-'(P»))2, 


If  we  assume  that  G(x)  = x for  0 < x < 1,  then  the  Hr,  (e)  is 


1 


given  by  (1.2.29)  with  a = u,  h(a  ) = 1,  n = ct  and  g(n  ) = 1. 


Hence 


np_  ( <z) 


(1.2.44)  ARE(R,,R.;6)  = lim  — 

I 5 €>o  "r/^) 


= B^d-g)  , 3(1-g)(2-g)^  ^ ^ 

^ 3(l-g) (2-g)^+a^ 


= 0.931,  a = 


_ 1 

2- 


= 0.675,  a = J. 


= 0.574,  a = 


(C)  Selection  procedure  for  distribution  -<  ordered  with  respect  to 

c 

Weibull  distribution 


Assumed  that  the  specified  distribution  G(x)  is  a Weibull 
distribution.  In  other  words,  G(x)  is  given  by 

■ax'' 


G(x)  = I 


1 -e 


0 


for  X 4 0 
for  X < 0 


jlA 
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where  > •»  0 and  attention  is  restricted  to  a ^ 1 which  is  assumed 
known. 

In  this  case,  we  use  the  statistic  Tt  which  is  defined  by 

r 


(1.2.45) 


X“ 


j = l 


1 ;j  .n  ;r,n’ 


i = 


(as  before,  X...  denote  the  j-th  order  statistic  from  F.,  i=l,...,k). 

I > J » n 1 

Since  G(x)  is  convex  with  respect  to  the  exponential  distribution  if 
a ^ 1 and  since  the  convex  ordering  is  transitive,  the  family  of 
distributions  which  are  convex  with  respect  to  Weibull  (a  ^ 1)  will 
have  IFR  distribution.  Thus  our  interest  here  is  in  a special  subclass 
of  IFR  distributions. 

The  rule  for  selecting  the  population  which  is  associated  with  F^j^j 
is  as  follows, 

Rg;  Select  population  if  and  only  if 

(1 .2,46)  Tt  > Cc  max  Tt 

' ■ ® l<j<k  J 


where  Cg  (0  < Cg  _<  1 ) is  determined  so  as  to  satisfy  the  basic 
probability  requirement. 

Using  the  fact  that  if  F < G and  F(0)  = G(0)  = 0 

c 


then 


where  F is  the  c.d.f. 

a 


c.d.f.  of  Y”'  and  G(y) 


F G for  a > 1 , 
oc  c « 

of  X“,  F(x)  is  the  c.d.f. 
is  the  c.d.f.  of  Y. 


of  X,  G is  the 

a 


Also, 


n^  St  order  statistic  from  G*(x)  = 1-e"^^ 


X > 0 where  X,  „ <...<  X„  „ are  order  statistics  from  F.  In  a 


■p^ 

I 

manner  similar  to  the  proof  of  Theorem  1.2.1,  one  can  prove  the 
following  theorem. 
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i 


1 


I 

i 

I 


Theorem  1.2.3.  IfF-(x)2^Fj-|^j(x)  and  i = 1 

Fj-^j(O)  = 0.  Fj-^j  <G,  G(x)  = l-e''""  . x > 0. 

r- 1 

is  known.  If  n max(r,  .j },  then 

'■^6 


. . . ,k, 

A > 0 and  u(^  1 ) 


(1.2.47)  inf  P[CS|Rj  = / G^'''(J-)dG  (x) 

° 0 6 

2 

where  Gj(x)  is  the  c.d.f.  of  a >:  random  variable  with  2r  d.f. 


( D ) Selection  with  respect  to  the  means  for  Gamma  populations 
Let  be  k populations  with  densities 


. i -.X 

r . ( X ) - -:rf \"X  e ,X  '0,i:>:’0,oi->l»i”l»...»k. 

1 1 1 '*  ' 

Let  R.j(x)  be  the  distribution  function  of  i = l,...,k.  We 

are  given  a sample  of  size  n from  each  n^.  Let  T*  be  total  life 

statistic  until  r-th  failure  from  -n.!. 

Let  aj-^j  1---1  be  the  ordered  values  of  a.j's.  We  are 

interested  in  selecting  the  population  with  the  largest  value 

^i 

(unknown).  Since  the  mean  of  is  — , selection  of  the 
population  with  largest  mean  is  equivalent  to  selecting  the 
population  with  largest  value,  ^be  subset  selection  rule 

based  on  T^  is: 


I 


R^:  Select  population 


if  and  only  if 
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(1.2.48)  Tf  > c,  max  T*. 

where  c-j  (0  < < 1)  is  chosen  to  satisfy 

PLCSlR^]  ^ P*. 

Since  the  rule  is  scale  invariant,  we  can  assume  6=1. 

Case  1:  All  a.  are  unknown  and  1 . Let  = {a  = (a^ ,. . . .aj^) : a^.  ^ 1 Vi} 

R^  ■<G(x)  = 1-e'^,  X - 0,  i = l,...,k. 
c 

In  this  case,  by  Corollary  1.2.2,  we  have  the  following  result. 

If  n > maxir,  f~},  then  inf  P[CS|R7]  = (^)dH(x) , 

' 0 ^ 

2 

where  H(x)  is  the  c.d.f.  of  a x r.v.  with  2r  d.f. 

Case  2:  a.  are  unknown  but  assume  1 i = l,..,k  and  A is 

known . 

Let  R^(x)  be  the  c.d.f.  of  X with  density  function 

r^(x)  = x'^'^e’’'^,  X > 0,  3 > 0. 

2 

Let  h(x)  be  the  density  function  of  a x r.v.  with  2r  d.f. 

Theorem  1.2.4. 

(1.2.49)  PCCSIRJ  i/  h''-'(^x)  e-^  dx, 

where  y = (1  . 

A 

* 

Proof.  PLCSIR^]  = P[T(^)  > C7^max^ 

where  is  associated  with  = l,...,k. 
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1 

I 


Since  R^(x)  £ R^-(x)  ^ G(x)  = by  Lemma  1.2.6,  then 

(1.2.50)  PCCSIR,]  > P[T**>  C-,  max  T**J 

where  Tj^*is  the  total  life  statistic  until  r-th  failure  from  G(x) 

and  T|*(j  = l,...,k-l)  is  the  total  life  statistic  until  r-th 

failure  from  R^(x). 

Since  a > 1 then  R.  •<  G. 

^ c 

Let  cf  (x)  = G'^R^(x). 

(1.2.51)  P[T**>  C7  TJ*,j  = l,..,k-l]  = P[cp(i  TjJ*)>  cp(^  Tj*)  j = l,...,k-l] 

c^  (n-r+l)c^ 

By  Lemma  1.2.4  (a)  with  a^  =...=  a^_^  = — , 

a^  = 0 for  i > r+1 


and  cf(x)“^  Y where  X is  a r.v.  with  distribution  R^ 
Visa  r.v.  with  distribution  G, 


we  have 

c c 

(1.2.52)  P[cf(l  T**)>  cf</TJ*).j=l,...,k-l]  > P[cK^  T**)i  ^ Y., 

j = 1 , . . . ,k-l] 

2 

where  Yj  (j  = l,...,k-l)  is  a r.v.  with  x with  2r  d.f.  From 
(1.2.50),  (1.2.51)  and  (1.2.52),  we  have 


PLCSlR,]  > / x)dB(x), 

' 0 n 

where  B(x)  = P[cf(-  TM);<  x]. 

n K 

Since  q(x)  = ( 1 -R^(x) ) , then  ^(x)  = pj(1-e  ^). 
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B(x)  = P[T**<  n cf‘^(x)]  = H(2n  cp“^(x)) 
= H[2n  R'^d-e-")]. 


^ = h[2ny]  ^ 


= h(2ny)2n 


where  y = (1-e  . 


So  that. 


1 x)dB(x)  = / X)  e-«dx. 

This  completes  the  proof. 

Let  S denote  the  size  of  the  selected  subset.  The  expected 
value  of  S when  Ry  is  used  is  given  by 

k 

(1.2.53)  E(S|RJ  = y P[TT  > c,  max  T^]. 

i = l lUlk  ^ 

Let  n’  = {a  - (a.|  , . . . ,a|^ ) : 1 £ £4,  i = l,...,k}. 

For  a € fl',  since  R^(x)  £ R.j(x)  £ G(x)  = 1-e  then 

k 

ECSIR,)  < I P[T**>  C-,  max  T^*] 

^ " i = l ' ' 2£j£k  ^ 

where  T|*is  the  total  life  statistic  until  r-th  failure  from  R^(x) 
and  TT*(j  = 2,...,k)  is  the  total  life  statistic  until  r-th 
fail ure  from  G(x) . 

Hence  E(S|R,)  < k P[T?*^  c,  max  T**J 
' ' ' 2<j<k  J 

so  that, 
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(1.2.54)  sup  E(SiRj  = k /R'^‘^(7:^)dS(x) 

I I * ' ^ “J 

where  R(x)  (S(x))  is  the  c.d.f.  of  total  life  statistic  until 
r-th  failure  from  G(x)  (R^(x)). 


1.3.  Selecting  a best  population  - using  indifference  zone  approach. 


Let  be  k populations.  The  random  variable  associ- 

ated with  TT^  has  an  absolutely  continuous  distribution  F^.  We  assume  there 
exists  a such  that  i = l,...k-l  and 

6 (0  < 6 < 1)  is  specified.  Let 

(1.3.1)  = {F  = (F^ ,. . . ,F|^) : aajsuch  that 


^i(x)  1 Fj(J)  V i M)  . 

The  correct  selection  is  the  choice  of  any  population  which  is 
associated  with  We  propose  the  selection  rule  Rg:  Select 

population  if  and  only  if 

(1.3.2)  T.  = max  T-  where  T.  is  defined  as  in  Section  1.2. 
^ llj<k  ^ ^ 

We  want  the  P[CS|Rg]  P*,  for  all  F fc  Q , where  P*  (^  < P*  < 1 ) 
is  specified. 


Theorem  1.3.1.  If  F^. , G € y,  i = l,...,k,  P^|^j(0) 
and  G”^  (0)  £ 0.  If  a^  ^0,  j = 1 , . . . ,r, 

gG'^  (0)  £ 1 and  a^  ^ 6, 


0, 


< G 
c 


then 
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(1.3.3)  P[CS|Rg]  > / Gif’(|)dG.p(x) 

where  Gj(x)  is  the  c.d.f.  of  T. 


Proof.  P[CS|Rg]=  > max 

Since  Fj-.j((Sx)  > Fj-(^j(x),  i = l,...,k-l,  it  follows  from  Lemma  1.2.6  that 


P[CS|Rg]  = P[T^^^  > 64“^^  V j f k] 


P[T(^)  > 6 Tj  vj?«k] 


where  are  i.i.d.  with  c.d.f.  Wj^(x).  Using 

the  same  argument  as  in  Theorem  1.2.1,  we  have  our  theorem. 


Remark  1.3.1.  inf  P[CS|R]  = / G^'^  (|)dG.j.(x)  if  G'^(O)  = 0. 

G'^(O) 


For  given  k,  6,  P*  and  G(x),  we  can  possibly  find  the  values  of 
the  pair  (n,r),  (n  ^ r)  which  satisfy 

(1.3.4)  a^  > r and  / g!J’^  (|)dG.j.(x)  > P*. 

G‘^0) 

If  G(x)  = X for  0 < X < 1 , we  can  always  find  the  values  of  the 
pair  (n,r),  (n  ^ r)  which  satisfy 

']  Glf’^f)dG,(x)  > P*. 

0 

In  this  case,  Gy(x)  is  the  c.d.f.  of  r-th  order  statistic  based  on  a 
sample  of  size  n from  G. 

If  G(x)  = 1-e'^  for  x > 0 (see  Patel  [58]),  we  can  find 
the  smallest  integer  r which  satisfies 


j 


F 
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o . _ 

[ G-p  (j)dG-p(x)  > P*  where  G-p(x)  is  the  c.d.f.  of  a random  variable 
° 1 

with  2r  d.f.  Since  ^ > 1 iff  n > ^ , the  minimum  n 

r- 1 

satisfies  n > max{r,  t— r-}- 


1.4.  Selection  procedure  for  distribution  -c  ordered  with  respect  to  G. 

★ 

Let  t(>  0)  be  a given  number.  Let  N^.  (t)  be  the  number  of  failures 
in  time  t among  the  n units  on  life  test  from  u^.  which  has  a 
continuous  distribution  , i = l,...,k.  We  assume  that 


(1.4.1) 


N^.(t)  > 1 , i = 


and  there  exists  a such  that  F^.(x)  ^ 

i = l,...,k.  The  correct  selection  is  the  choice  of  any 
population  which  is  associated  with  Let  T.j  be  the  total  life 

statistic  until  N.(t)-th  failure  from  population  . Let  ^ be 
associated  with  F|-^.  j.  We  propose  the  rule 

(1.4.2)  Rq:  Select  tt.  » T.  > Cq  max  T. 

lUlk 

where  Cg(0  < Cg  £ 1 ) is  chosen  so  that  P[CS|Rg]  £ P*. 


Theorem 


1-4.1.  If  F^.  (x)  ^ F|-|^j(x)  y X,  i = 1 ,. . . ,k  and  F|-|^j  < G, 


(1.4.3)  P[CSiRq]  > / A'f-''(MdA.(x) 

^ 0 9*^ 

where  A-j  (x)  (A2(x) ) is  the  c.d.f.  of  total  life  statistic  until  n-th 
(first)  failure  from  G. 
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Proof.  From  (1.4.1),  we  have 

P[CS|Rg]  = > Cg  max 

■ --  ^9  j " l,...,k-l] 


where  — - is  the  mean  of  a sample  of  size  n from  i = l>- 

and  — is  the  first  order  statistic  from  P[-|<.-]- 


.,k-l 


(1.4.4) 


Tr  = I X.  for  j = l,...,k-l. 


where  Xj^,...,Xj^  are  i.i.d.  from  j = l,...,k-l. 

sine  1 'I'j-  Hence 

(1.4.5)  P[T*  > CqTJ,  jfk]  > P[T*  > CgTJ*,  j^k]. 

Let  q^x)  = G‘V|-^j(x). 

T* 

(1.4.6)  P[T*  > CgTj*,  jikj  = P[c?(~)  i Tj*),  j^] 
Since  c,  is  starshaped,  then 


(1.4.7) 


T*  T* 

/]_  . 1_  {A\ 

'^^Cg  n ^ - Cg  "^n  '■ 


By  Lemma  1.2.5  (a)  and  (1.4.4), 


(1.4.8) 


4 Tf ) ' I .<X^, 


'k 

Note  that  cf(— -)  Y|^  where  Y^,  is  the  first  order  statistic  of 

size  n from  G and  has  distribution  G,  j = l,...,k-l. 


?,  = 1 , . . . ,n. 
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Let  Y = I cp(X.  J,  j = 

From  (1.4.7)  and  (1.4.8),  the  right-hand  side  of  (1.4.6)  is  greater  than 
or  equal  to 

4,  ’k  J'l.-'-.k-i:  ' i''j.  i=' w] 

rx) 

= { A!f‘^(^)dA.(x). 

0 ^9 

The  ..  .of  is  complete. 


1.5.  Selection  with  respect  to  the  means  for  distributions  s-ordered 
with  respect  to  a specified  distribution  G. 

Let  the  random  variable  X^.  associated  with  have  an  absolutely 
continuous  distribution  F..  Assume  that  F.(x)  = F(|— ),  i = 1,...,k, 

7 7 e. 

F < G and  F and  G are  symmetric  about  0.  As  before,  let 
s 

0 ■■  Opj  be  ordered  values  of  unknown  e^.  (i  = 1,...,k). 

A best  population  is  defined  as  the  one  having  We  are  given 

a sample  of  size  n from  each  of  the  k populations.  Assume  that  n 
is  large.  Let 

('■5-D  T,  . X,_,,yx,.^.^,^„.... 

where  X,  X.  “ i *i  ;t,-H  ,n  ' 'i  5 ;n,n 

order  statistics  from  F^(x),  and  r',r  are  given  integers  such 
that  1 < r'  < p,. , n ' r > , i = 1 , . . . ,k. 

Similarly,  let 


I 

'( 
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(1.5.2)  T = -r'Y  , -Y  -...-Y  +Y  +...+( n-r+1 ) Y 

^ ' r ,n  r +1 ,n  . ,n  r+1 ,n  ' ' r,n 

where  Y,  Y 0 ' Y Y„  „ are  order  statistics 

1 ,n  “ - 1, ,n  /.+1  ,n  - - n,n 

from  G.  Then,  for  selecting  a subset  containing  we  propose 

the  following  selection  rule: 

Select  population  if  and  only  if 

(1.5.3)  T.  yc,^  max  T . 

where  c-jg  (0  ■'  c.|g  1)  is  determined  as  to  satisfy  the  probability 
requirement,  inf  P[CS|Rig]  ■ P*,  where  CS  stands  for  a correct 
selection,  i.e.,  the  selection  of  any  subset  which  contains  the 
population  with  and  ,2={  (o^ ,. . . ,0|^) : O^.  > 0,  i = l,...,k  . 

Let  j be  associated  with  i = l,...,k.  Now,  we  give 

a theorem  pertaining  to  the  computation  of  the  minimum  probability 
of  a correct  selection. 


Theorem  1.5.1.  If  F.  (x)  = F(— ) vxandi  = l,...,k,  F<Gand 

^ '■’i  s 

F and  G are  symmetric  about  0. 

If  = t',  lim^  = t and  c^g  • min'f  , 1-ti, 

n-MX)  n 


then , 


(1.5.4)  lim  inf  PLCSjR,.]  = / )dH(x) 

n-K.  s,  '^0  ^10 


where  H(x)  is  the  c.d.f.  of  T defined  in  1.5.2. 


Proof.  P[CSiR^g]  = P[T(^j  ^ ^io'^(j)’  J = 1 • >>^-1 1 

whGrGT/*\  “ **r  X/*\  i i.t  f a \ j.i 

(j)  i^)r\n  (i)r'  + l ,n  (i  (i ) w • ^ + 1 ,n 


(i)’"  ^'^'(i) 

+ (n-r+l  )X 


(i  )r,n 
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and  X, 


< 0 < X,. 


<...<  X, 


order  statistics  from 

It  is  easy  to  see  that 


are 


(1.5.5)  inf  P[CS|R^q]  = P[T-  > c^qH,  j = l,...,k-l] 

where  T'.  = -r'Y.  , -...-Y.  „ +Y.  ,,  +. . . + (n-r  + l )Y . 

1 T,r',n  iU^.j,n  iu^.^  + l,n  i;r,n 

3nd  Y-  1 „ <...<  Y.  . „ " 0 < Y-  „ „ <...<  Y.  „ „ are  order 

statistics  from  F(x) . 

Let  ■j.(x)  = G’ V(x) . 

n.5.6)  P[T'  , c,„Tj]  . PU(jL  rp  > j = I,,..,k-l]. 

By  Corollary  2.9  of  Lawrence  [48]  with  ig  = Jq  = z^.^+1, 

and 


then 

C-5.7)  .(J-t;)  : — .....  J 


Also  by  Corollary  2.10  of  Lawrence  [48]  with  i„  = r'-l,  = n+1 , 


0 

0 •_  j 

-r' 

j = r' 

-1 

ir  * ' 1 

' J 

_1_ 

"^10 

'■(k)  ^ 

n-r+1 

"^10 

j = n 

V » 

< .1  and 

nc^O  - nc^Q 


"(k) 
j r 


J > r 


k;j,k,,n' 


— ^(Y,  1 

"‘^10 


)+...+ IL:ni^(Y^  ). 

' ncm  k;r,n' 


is  the  i-th  order  statistic  from  G(x),  j = l,...,k. 


j'.'  = -r'V.  , -...-V.  +V.  1 +. 

1 T,r',n  ''’'(i)’'^  TU(,-)+l>n 

From  (1.5.6),  (1.5.7)  and  (1.5.8),  we  obtain 
n-^i  1 c^qTJ,  jfk]  > P[T|;  ic^plj,  j = l,...,k-l]. 
Since  G is  symmetric  about  0,  as  n becomes  large, 


3,T,  j = l,...,k. 


Thus 


P[CS|Rio]  1 / H''‘''(|-^)dH(x)  if  c^Q  < ■ 

This  proves  the  theorem. 


-r'  n-r+1 


Corollary  1 .5.1 . If  F-(x)  = F(^)  v x and  i = 1 k,  F <G,  F is 

"i  s 

symmetric  about  0 and  g(x)  = = -je'^^L-®<x<«=. 

r ' 1 r 1 

Iflim  — =t'  <2»  lim-=t>^  and  c.,^  ;<  min{t‘ , 1-t}, 

n-KX) 


n-H» 


then,  as  n is  large. 


10 
1-c 


(1.5.9)  inf  P[CS|R.  J «./  +(-r^)/F^PTT  ]d$(x) 

n — ‘^lO 

Proof.  Considering 

where  Y,  <...<  Y.  „ < 0 < Y„.,  „ <...<  Y_  . are  order 
i,n—  — t,n  ?.+  l,n  n,n 

statistics  from  G(x) . 

2 

Since  2 T*  is  distributed  as  a x >^-v.  with  2(r-r'+l)  d.f. 
(see  Lawrence  [47])  and  T*  has  same  distribution  as  T when  n is 
large,  then  2 T is  distributed  as  with  2(r-r'+l)  d.f. 

By  Theorem  1.5.1, 

P[CS|R^q]  > P[V|^  > C^g  max  Vj] 

^ 1 

2 

where  V.|,...,V|^  are  i.i.d.  from  x with  2(r-r'+  1)  d.f. 

Since  - N(2(r-r'+l),  4(r-r'+l)),  then 

P[V.  ^ c max  V .]  - / (^^ — + ^ /r-r  ' + 1 )d<t(x) 

I j k J -■  ^10 

This  proves  the  corollary. 


■10 


w 
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Assume  ■ -1  ^ 5 0 ^k+1 -**•-  n order  statistics 

from  distribution  G where  ^ G(x)  = g(x)  = ^ e*  ^ ^ < x < =“. 

Define 


i-1 


(1.5.10)  V, 


. ifi 


n 'n' 


, i = 1 


,n, 


Where  (i)  = gG'^J^)  (X^  ^ J . 

Now  we  want  to  get  the  limiting  distribution  of  (as  n is  large). 

Let  = /n[l  - ^].  To  obtain  the  asymptotic  distribution  of 

we  use  the  following  result.  Let  X,  <...<  X „ be  the  order  statistics 
^ 1 ,n  n ,n 

from  G and  S„  = ^ L(1)X,_„. 


/ = a^(G)  = 2 / / L(G(s))L(G(t))G(s)[l-G(t)]dsdt. 
s<t 

where  L satisfies  the  condition  (ii)  of  following  lemma. 


Lemma  1.5.1.  (Moore  [53]) 

2 

If  a < '»  and 

1 T 

(i)  E|X|  = /|G’l(u)|du  < - 
0 

(ii)  L is  continuous  on  [0,1]  except  for  jump  discontinuities 
at  a.|,...,a|^^,  and  L'  is  continuous  and  of  bounded  variation  on 

[0,1]  ~ {ai,...,a|y|}, 

then 


(1.5.12)  f/n  [S^  - / xL(G(x))dG(x)]}  - N(0,a^). 


Let 
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Then  using  an  argument  similar  to  the  one  in  Barlow  and  Doksum  [3] 
1 " 


^.n 


where  L(i)  = g G“'(i^)-n(i  - ^)[g  G-^i)-g 

1 

g ’g(x)  = 


Since 


for  0 < X £ ^ 


1-x  for  ^ < X £ 1 


then 


'■(^1  = 1 


2(1)  - } for  1 < ^ 


I ■ fo'-  M 


Thus  we  can  consider 
(1.5.12)  S 


J,  "-(ijXl.n 
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(1.5.13)  a^(G)  = 2 / / L(G(s))L(G(t))G(s)(l-G(t))dsdt. 

s<t 


2(G)  = 2 sds]  --  4?-^-  . 0-t)dt 


1 JUO J - 1 

g G '(s)  gG'(t) 


= 2 


i t 


I [/  sds]  (l-t)dt 

6 0 g G‘‘(s)  g G"'(t) 

+ 2 / [/  — sds  + / sds]  dt 

X 0 gG"'(s)  X gG"'(s)  g G"  (t) 


2 . J + 2 • ^ ^ 

^ 48  ^ ^ ■ 4 • 


By  Theorem  2.1  of  Barlow  and  Van  Zwet  [11],  H~^(l)  G(X  ) a.s. 

n n ^ D 

Since  G(X^  ^)  1,  therefore  H^^(l)  >■  1 a.s. 

/nS 

Now  T = — T , hence  by  Slutsky's  theorem, 

" h;;Vi) 

we  have  the  following  result. 


Theorem  1.5.2. 

(1-5.14)  ^ - 1))  vN(0,  \) 

1.6.  Selecting  a subset  which  contains  all  populations  better  than 
a control . 

In  addition  to  TTp...,-n|^,  we  have  a control  population  ttq.  Let 
have  a continuous  distribution  , i = 0,1, ...,k.  We  are  given 
a sample  size  n from  each  of  the  (k+1)  populations  tt^,  i = 0,1,..., k. 
Let  be  the  unique  r, -quantile,  i = 0,l,...,k.  A population  n. 

(i  = l,2,...,k)  is  defined  as  better  than  the  control  ttq  if  ■ f,^g. 
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i 

i 


i 


1 


We  want  to  select  a subset  such  that  the  probability  is  at  least 
p*  (1  < p*  < 1 ) that  all  populations  (i  = l,...,k)  for  which 
^ ^ will  be  included  in  the  subset.  We  regard  any  such 

selection  a correct  selection  (CS).  We  assume  that  for  a fixed 
i , 1 = 1 , . . . ,k,  either  (x)  1 Fq(x)  or  F^.  (x)  1 Fq(x)  V x > 0. 

Thus  Ti^  is  better  than  iff  F^(x)  5 V x > 0.  Here 

tiq  is  not  known  and  we  assume  Fq  < 6.  Let  T.  be  the  r-th  order 
statistic  from  F^  where  r ;<  (n+l)a  < r+1  and  let  H.  be  the  c.d.f. 
of  the  r-th  order  statistic  from  F^. , i = 0,1, ...,k.  We  propose  the 
rule  R.|i  as  follows, 

Select  tt^-  if  and  only  if 

(1.6.1)  T.  > c^^Tq 

where  c^^  (0  < c^^  < 1)  is  determined  to  satisfy  the  basic  probabil- 
ity requirement.  Let  t be  the  number  of  tt^'s  better  than  tiq  and  let 
A be  the  index  set  for  those  , t is  unknown.  Let  G^(x)  be  the 
c.d.f.  of  the  r-th  order  statistic  of  sample  size  r from  G.  Let  ii 
be  the  set  of  (k+l)-tuples  (Fq,F^.  . . ,Fj^). 

Theorem  1.6.1.  If  F„  -<  G and  for  a fixed  i,  i = 1 k,  either 

^ - u * 

F.(x)  > Fq(x)  or  F.(x)  < Fq(x),  x > 0,  then 

(1.6.2)  inf  P[CS|R.|^]  = fc(c^^x)dG^(x)  where  5^(x)=l-6^(x) . 

a 0 

Proof.  P[CS|R^^]  = P[T.  >c.|^  Tq  i e a] 


n [l-H.(c^^x)]dHQ(x) 


I 
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Since  F^(x)  ^ ^ H^.(x)  ^ Hq(x). 

Thus 

«x> 

P[CS|R,,]>/  TT  [l-H.(c,,x)]dH.(x) 

“ 0 1 a ^ ^ 


= / [I-Hq(c^ix)]^  (IHqCx)  [l-HQ(c^^x)]'^dHQ(x) 

= P[  min  Z.  > c ^ Zg] 

1 <1  <k 


where  Zg,  Zi,...,Z|^  are  i.i.d.  with  c.d.f.  Hg(x).  Since  Fq  < G 
(see  Barlow  and  Gupta  [5]),  then 


P[  min  Z.  ^ c,^  Zg]  > /[6^(c,  ,x)]*^  dG  (x). 
1^i_<k  0 

This  completes  the  proof. 


1.7.  Quantile  selection  rule  for  a restricted  family  of  distributions 

We  assume  that  each  F^(G)  has  a unique  a-quantile  ^nd 

F.(0)  = 0 = G(0),  i = 1,2,.. .,k. 

Let  Cm  <...<  CruT  be  ordered  values  of  c,-  . i = 1 > • • • .k.  A 
L I Ja  — — [k Ju  la 

best  population  is  defined  as  the  one  having  Let  T^. 

denote  the  r-th  order  statistic  from  F^. , where  r (n+l)a  < r+1 . 

The  rule  we  propose  is 

R-J2:  Select  if  and  only  if 

(1.7.1)  T.  > max  T . - c,, 

’ l<j<k  J 

where  c-|2  (0  < c-|2  1 1 ) is  determined  so  as  to  satisfy  the  probabil- 
ity requirement  P[CS|R^2l  1 P*- 


L 


As  before  let  Q be  the  set  of  all  k-tuples  (F-|,.,,F|^).  Also 

if  F has  a unique  a-quantile  ^ , we  define  F •<  G if  and  only 

“ S* 

if  F(0)  = 0 = G(0)  and  G”V(x)  is  concave-convex  about  the  C , 

a 

on  the  support  of  F. 

For  example,  if  G(x)  = 1-e"^,  then  F < G is  equivalent  to 

S* 

saying  that  F(x)  is  DFR  for  x < ^ and  F(x)  is  IFR  for  x > c . 

— a a 

on  the  support  of  F.  We  are  interested  in  this  "turning  point" 

f,  . Let  H .(x)  be  the  c.d.f.  of  the  r-th  order  statistic  from 
a r , 1 

F|-^j  and  let  G^(x)  be  the  c.d.f.  of  the  r-th  order  statistic  from 

G. 


Theorem  1.7.1.  If  ^ ^ i 0 i = l,...,k, 

""[k]  ^ i G(y+nJ,  for  y > 0,  then 


(1.7.2) 


Inf  P[CS|R^2^  = / Gjl*''(x+c^2)dG^(x) 


where  G (x)  is  the  cdf  of  r-th  order  statistic  from  G. 
r 


Proof. 

(1.7.3) 


k-1 


P[CSiR,2l  = ^ 


,k-l 


12 


where  X^^...,X^  are  i.i.d.  r.v.'s  with  distribution  |^. 


Let  c,  (x)  = G*'Fj-|^j(x)  = G^'^r  k^^^'  concave- 

convex  about  ff'(x)  is  decreasing  in  x for  0 < x < c 


[k]  ,01 
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and  q:'(x)  is  increasing  in  x for  x > r 


[k]a- 


So  that  cf'(x)  > tf' V X > 0.  Since  1 G(y+n^) 

for  , . 0,  then  F|-^^(x)  i for  x , 


Hence 


- "o  , 


x-f 


[k]a 


Let  A(x)  - G - G = G - n^,  then 


"‘"^[klc 


Let  X 


Si  nee 


'■[k]n-  ''■(f-[k]«>  i '• 


dx 


^rirT(x) 


- A(x)  = A'(x)  = M 


g[G‘^F[^](x)] 


then 


Hence  cf'  (x)  ^ 1 v x > 0. 


Si  nee 


c(<max  X^^.)  - cf(x^^, 

max  X . - X ! 
r,i  r,k 


cj'(t)  for  some  t such  that 


X„  1,  1 t j<  max  X . and  f4{X  ,• ) = V . has  distribution  G , 
r,K  r,i  r,i  r,i  r 

then 


^r  i ■ ^r  k - Vi'  V k' 
l<i<k  l<i<k 


From  (1.7.3),  we  obtain 


P[CS|R,j]  > P[^max^  < c,^] 


= / G^'^(x+c,2)‘1G^(x). 


This  completes  the  proof. 


1.8.  Estimation  of  ordered  parameters 


We  assume  that  (x)  = F(^)  vx  and  i = l,...,k,  where  F 

and  0^  are  unknown  and  • 0,  i = l,...,k.  Let 

®[1]  1- • -1  6[|<;]  ordered  values  of  6^'s.  We  assume  that 

each  F^(F)  has  a unique  a-quantile  i = l,...,k. 


Let 


‘^[1] 


< . . . < 


be  ordered  values  of  C^'s.  We  are  given  a 


sample  of  size  n from  each  of  the  k populations.  Let  T^.  be 
the  total  life  statistic  until  r-th  failure  from  F^. , i = l,...,k 
T*(T)  be  the  total  life  statistic  until  r-th  failure  from  F(G) 
and  let  Tt  be  the  r-th 
where  r < (n+1 )a  < r+1 


and  let  Tt  be  the  r-th  order  statistic  from  F^. , i = l,...,k. 


Let  1- ■ -1  be  order  values  of  T^. 's 

and  let  be  order  values  of  Tf's.  Let  and 

be  associated  with  i " l,...,k. 

(A)  Estimation  of  Oj-^j  based  on 

Theorem  1.8.1.  If  F < G,  F(0)  = G(0)  = 0 and  if 

c 

00  00 

9n  = f xdF(x)  = / xdG(x)  where  e,-,  is  known,  then 
^0  0 

(1.8.1)  ^^®[i]  — ^[i]^  — ^ ^ 

_I 

where  a = W*^(1  - and  W(x)  is  the  c.d.f.  of  T from  G. 


- A ^[i]^  " ^ ®[i]^- 


Proof. 
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I 


Since  i t]  is  a non-increasing  function  of  (1  ^ k) 

(see  Chen  and  Dudewicz  [20  ]).  then 


= 1 A]. 

By  Theorem  4.6  of  Barlow  and  Proschan  [8], 

p‘^-^^^[T*  > a]  > [T  ^ a] 

1 

= n-w(A))^-^"^  = n-(i-,'^),k-iei  ^ 

2 

Let  x,^(2r)  denote  the  100a  percent  point  of  a chi-square  distribution 
with  2r  d.f.  Again,  by  Theorems  3.3  and  2.4  of  Barlow  and  Proschan 
[9],  we  have  the  following  theorems. 

a> 

Theorem  1.8.2.  If  F is  IFRA  and  / xdF(x)  = where  0^  is  known. 


(1.8.2)  ^l''’[7]^  -'i 

2 

{ -'O^l-x’^^r) 


y^_  , (2r)  " 2(n-r+l ) 


xf_p.(2r)  2(n-r+l) 


where 


k-i+1 


Theorem  1.8.3.  If  F is  IFR  and  / xdF(x)  = 9^.  where  60  is  known. 


then 
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(1.8.3)  ^ c^Tj-^.jJ  > t>2  where  c^ 


2nk 

OoXgn(2r) 


1 

and  f-'"  = ^2  ■ 

Let  ByP.2  (O  ■'  B-|,L2  1 ) be  such  that  B-|  + B2~l  = y where  y is 

fixed  and  Q < y < }.  By  Theorems  1.8.2,  and  1.8.3  and  the  fact  that 
P(AB)  P(A)+P(B)-1  with  A = ^^[ij  — snd  B = — ^2"^[i]^’ 

we  have  the  following  theorem  for  two-sided  confidence  interval 
of 


Theorem  1.8.4.  If  F is  IFR  and  f xdF(x)  = 6n  where  6,^,  is  known, 

0 u u 

then 


(1.8.4)  P[CjTfjj  < e^.j  <_c,T[^j]  V 
where  c.j , C2  are  defined  as  above. 

(B)  Estimation  of  based  on  Tj-^j 
Theorem  1.8.5.  If  F ^ G,  then 


\ 


(1.8.6) 


[ 


where  A.j 


if  (n-r+1 ) 1 

if  A-]  (n-r+1 ) £ 1 


G'h.,) 


and 


W(x)  is  the  c.d.f.  of  T from  G. 
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Proof.  iC,T|-,j] 


,k-i+1 


Theorem  1.8.1) 


f^k-i+lfT  ^ '"1 

fTi) --  c7  ■ "[i]  ’ 

" ?*"-'♦' [T*  ^ J-  ,-] 

^1 

= > ^-F-^c)]  = p'^‘'^^[F(c,T*)  > a]. 

Cl  I 

In  a manner  similar  to  the  proof  of  Theorem  3.1  of  Barlow  and 

Proschan  [9],  we  can  show  that 

1 


P[F(c^T*)  > > B, 


k-T+T 


Hence 


Looking  at  Theorem  2.1  of  Barlow  and  Proschan  [9],  we  can 
have  the  following  theorem. 


Theorem  1.8.6.  If  F G and  F(0)  = 0 = G(0),  then 

( 


(1.8.7)  > '^2  ^[i]^  - ‘■'2  ^2 


if  A^n 


i f A^n  1 


•1 


Oj  = -5_Lji  and 

u-’iT) 


W(x)  is  the  c.d.f.  of  T from  G. 

Let  P.|,  ' ''i  ’ ‘■'2  ' 1 "2'^  ^ ' where  , is 

fixed  and  0 ■ y ' 1.  By  Theorem  1.8.5  and  Theorem  1.8.6,  and  the  fact  that 

if  F(0)  = 0 = G ^(0)  = 0,  then  F ^G=^F<Gwe 


have  the  following  theorem  for  two-sided  confidence  interval  of 

^[i]* 
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Theorem  1 .8.7.  If  F < G and  F(0)  = 0 = G*^(0),  then 


(1.8.8)  1 1 


(c)  Estimation  of  Cr,-i  based  on  Tt 

nj  [ij 


Theorem  1.8.8.  If  F -<G,  then 

★ 


1 


(1.8.9)  (1)  P[^  < -M]  ^ B if  A,  = J ^ <1 

^1  ' G’’(a) 

1 

*T* 

T*  G • ( ■ ^ ) 

(2)  P['. r --|  ^ B if  Ao  ■ — T > 1 


■[i]  - A. 

Proof.  First,  we  want  to  prove  (1). 

pCH,] 


[i]  --"I'^ri]^ 

^ pk-i  + l ^ 1 Aif;j---j]  (see  Theorem  1.8.5) 


= fi-P[T*i)  ::  AiTj-.j]} 


k-i  + l 


k-i+l 


where  F.(x)  is  the  c.d.f.  'f  the  j-th  oi^Jer  statistic  from  F 

sJ 


Since  G"'F(x)  = g;)V-(x),  it  follows  that  F.(x)  = G-G"'f(x). 

J J J 


-1. 


Hence 


Fj(A^O  = GjG''F(A^r) 


; G.[a-,G’  F(0]  (since  F < G and  0 < A,  < 1 ) 
J ' * ' 

= Gj[A^G'^,)] 


= 1-e 


k-i  + l 


so  that 


Now,  we  want  to  prove  (2). 


vt.i’ 


> P’[T(i)  < A2^,|-ii] 

= f](a20 

= {GjG’V(A2:)}'' 

^ G.[ApG‘V(f,)]  t’  (since  F < G and  ‘.o  ^ 1 ) 

\_ 

= {Gj(A2G'\a))»'  = iGjGT^:- ')■*'•  = 

This  completes  the  proof. 
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CHAPTER  II 

INTERVAL  ESTIMATION  AND  SELECTION  PROCEDURES 
FOR  A SET  OF  GOOD  POPULATIONS 


2.1.  Introduction 


In  practice,  one  is  always  faced  with  the  problem  of  selecting 
the  better  ones  from  a group  of  populations.  Here,  we  consider 
such  kind  of  problems.  Suppose  that  we  are  given  k(k  ^ 2)  popula- 
tions -rr^ , . . . ,Ti|^  with  distributions  F(x,o.)  where  6.  lies  in  an 
interval  A on  the  real  line,  i = l,...,k.  The  quality  of  the  i-th 
population  is  characterized  by  the  real-valued  parameter  6.j.  The 
population  with  the  largest  0-va1ue  is  called  the  best  population. 

A population  is  considered  a good  one  if  its  quality  does  not  fall 
too  much  below  that  of  the  best  population.  If  d{e^.  ,9j)  is  a 
suitable  distance  measure  between  and  and  if  6|-|^j  = 
max  (0-1 ,. . . ,9|^) , population  r.  is 


(2.1.1) 


good 

, if 

d(r, 

[k]-''i 

bad. 

if 

d(o 

[k]’'i 

,i  ve 

constant. 

Let  X. 

i = l,...,k  be  mutually  independent  random  samples,  each  of  size  n, 
from  population  , i = l,...,k  respectively.  Let  T^.  " T(X^)  be 
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an  dpprotiridtp  p'.ti'njtor  of  . Let  q and  G (t,o.)  be  the 

1 1 n 1 

density  function  and  the  distribution  function  of  , respectively. 
The  results  of  tni-,  cruu.'ter  relate  to  two  cases, 


(i)  and 

(2.1.2) 

In)  G„;-,  n - t 0,  ...  0. 

If  the  distribution  function  F(x,n^)  of  belongs  to  a 
location  parameter  family,  then  F(x,t‘^)  = F(x-e^.).  The  distribution 
function  G (t,  .:)  of  i which  is  based  on  a random  sample  size  n is 

n I » 

not  necessary  of  the  form  G.(t-'‘^),  but  we  can  always  find  a 

statistic  T.  such  that  T.  has  the  distribution  function  of  the 
1 1 

form  G^(t--.).  For  example,  let  ii,,...,Tij^  be  k normal  populations 

p 

with  unknown  means  i'.] , . . . ,0j^  and  a common  known  variance  o '.  Let 
T^.  = be  the  sample  mean  of  based  on  a sample  of  size  n.  In 
this  case,  G^(x,h^)  (the  c.d.f.  of  X-)  is  of  the  form  G^(x,v  .j)  = 

G (x-h.),  where  G (x)  = f(—  x). 

ill  M O 

If  the  distribution  function  F(x,tm)  of  belongs  to  a scale 
parameter  family,  then  F(x,0^)  = F(^),  x • 0,  0^-  > 0.  As  before 
we  can  find  a statistic  such  that  has  the  distribution  func- 
tion of  the  form  G (|—) . The  results  of  this  chapter  deal  with 

n 

the  cases  when  G^(t,u^. ) is  one  of  the  two  forms  in  (2.1.2). 

In  Definition  (2.1.1)  we  take  d as  d^  or  d^,  respectively, 
whenever  o^.  is  a location  parameter  or  is  a scale  parameter  for 
the  density  of  T^. . dj^  and  d^  are  defined  as  d^(a,b)  = a-b, 
d^(a,b)  = 
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Estimation  of  based  on  (Tp...,Tj^)  has  been  considered  by 
several  authors.  Construction  of  two-sided  and  one-sided  confidence 


intervals  for  based  on 


(2.1.1) 


Trki  = max  T. 
l<i<k  ^ 


was  considered  by  Dudewicz  [24],  Dudewicz  and  Tong  [26],  Saxena  and 
Tong  [65]  and  others. 

For  problem  of  selecting  the  good  populations  from  a given 
collection  was  considered  by  Desu  [21],  Carroll,  Gupta  and  Huang 
[ 19]  and  others. 

We  consider  confidence  interval  I which  is  based  on 
In  this  chapter,  we  are  interested  in  finding  the  smallest  sample 
size  N such  that  the  probability  that  I contains  at  least  one  good 
population  is  at  least  P*  where  P*  is  a specified  number,  0 < P*  < 1 . 
The  probability  that  I contains  all  good  populations  and  the  proba- 
bility that  I excludes  all  the  bad  population  are  discussed. 

In  Section  2.2,  we  investigate  the  above  problems  for  location 
parameter  case.  The  infima  of  coverage  probabilities  are  obtained. 
Some  special  cases  are  discussed.  In  Section  2.3,  we  investigate 
the  above  problems  for  scale  parameter  case.  In  Section  4,  we 
discuss  the  above  problem  for  the  means  of  normal  populations  and 
for  the  scale  parameters  of  the  gamma  populations. 


2.2.  Results  for  the  location  parameter  case. 


In  this  section,  we  assume  that  is  a location  parameter 
for  T^ , i = l,...,k.  It  is  assumed  that  there  is  no  a priori 
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knowledge  about  the  o^'s.  For  arbitrary  fixed  d^ , d^  3 ^ 

confidence  interval  is  defined  as 


(2.2.1) 


M ^[k]  ^ ^2^ 


where  T 


[k] 


max  T,. 
1 •-  i < k 


For  the  location  parameter  case,  we  say  that  population  it^.  is 

(2.2.2)  good  if  2l 

bad  if  0^.  < 

where  eri  i = niax  e . . 

For  given  e = (0^...,^^^)  and  A > 0,  let  cx,j(l)  denote  the 
probability  that  I.|  contains  at  least  one  good  population,  let 
t.(2)  denote  the  probability  that  I.|  contains  all  good  populations 
and  let  a„(3)  denote  the  probability  that  I,  excludes  all  the  bad 


‘9 

populations. 


Let  us  denote  the  ordered  o-values  by  0|-.|j  5..  ..£  and 
let  ',1  be  the  parameter  space  which  is  the  collection  of  all  possible 
parameter  vectors  0 = («.| , . . . ,0|^) . For  given  a > 0,  let  m denote 
the  unknown  number  of  bad  populations  in  the  given  collection  of  k 
populations.  Clearly  we  have  0 • m < k-1 . 

Let 

(2.2.3)  = 10;  <...<  <...! 

Then 

k-1 

(2.2.4)  . 

.n-n  ’’1 

m-U 

We  need  the  following  lemma  before  we  prove  some  theorems. 
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Lemma  2.2.1.  Define  f(r)  = Ac'^-BD*^,  r = 0,1, ...,m  where  A 
If  either  (i)  1 < C £ D or  (ii)  1 > C > D > 0,  then 


B > 0. 


(2.2.4)  min  f(r)  = min[f (0) ,f (m) ]. 
0<r<m 


Proof.  Assume  1 < C ;<  D.  For  r = 0,l,...,m-l, 
f(r+l)-f(r)  = AC'^(C-1)-BD^{D-1) 

f(r+l)  £ f(r)  if  and  only  if  | ^ , r = 0,1 ,. . . ,m-l . 


Si  r,ce 
hence 


C ^ 1 A/Cxr+1 


D-1 


f (r+2)  < f ( r+1 ) . 


Similarly,  one  can  show  that  if  f(r+l)  > f(r)  then  f(r)  £ f(r-l). 
Hence  the  lemma  follows. 

If  1 > C £ D > 0,  then  by  a similar  argument  as  above,  the 
result  follows. 

This  completes  the  proof  of  the  lemma. 

Now  we  want  to  discuss  the  probability  that  I.|  contains  at 
least  one  good  population  and  the  infimum  of  this  probability  over 
We  need  the  following  lemma. 


Lemma  2.2.2.  If  the  family  of  density  functions  {g^(t-o):  o C 1 
has  a monotone  likelihood  ratio,  for  every  n £ 1,  then  for  arbitrary 
fixed  d.|  ,d^  satisfying  d.j+d2  > a. 


(2.2.5)  infa,,(1) 

' m 


min{G|I^(A+d^  )G|^'"'(d^  )-G^(A-d2)G|^’^  (-d^) , 


pk-m/ . 
^n 


-G|j''^(-d2) 
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Proof.  For  o c ;>  , 
m 


(2.2.6)  (,J1)  = P'^T|-|^j-d.|  ■ Oj-^. j < Tj-|^j+d2,  for  some  i,  i=m+1,. 

= l-P{0r.-i  (|;  I,  for  each  i = m+l,...,k} 


1-P{T[|^]+d2  --  ~ "^[k]' 

d^j+d^  > a) 

= P{T^I^]-di  < 0^^^}  - PfT^i^^+d^  < 


Let  t.  = 0[k]-^[i]^  i = l,...,k,  then  0 = 5^  < 5|^_^ 


'i  < < ...  < A 1 . 

m 1 


Hence 


(2.2.7) 


It  is  easy  to  see  that  a^(1)  is  nondecreasing  in  <5^+-].  where 
'm+2  ^ '■'m+1  '^m+l  ' ‘^m+2’ 

j=fm+2 


jfm+2 


In  this  case,  also  we  can  see  that  ji^{1)  is  nondecreasing  in  6 


where  6 ,,  < 6 


m+3  ^m+2  " ^et 


ni+2  ^ '^'tn+3'  Repeating  the  process  with 


each  6-,  j = m+l,...,k-l,  thus  we  have 

vJ 


(2.2.8)  Q inf  a^(l)  = :i  G^(.  +d^)C;;-"’(d^)-  n G„(^ ,--d,)Gr^-d,) 


where  = f ( , • • • .^^-1  ^ ^ ° ^ '^k-1  ■•••• 
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For  i = 1 ,2, . . . ,m. 

Hr  ' 5(a,+d,HA-B 

k-m  ^ k m 

where  A = G„  '^(d,)  n G f6.+d, ) and  B = G„  "’(-d,)  n G (5,-d,). 
n ij-inji  n 


g(6.-d2) 

Since  A > B and  ^ is  non-increasing  in  6^  > A,  i = 

then  Q is  greater  than  or  equal  to  the  right  hand  side  of  equation 

(2.2.8)  with  either  5^  = a or  6.j  = «>,  i = l,..,m. 

Tnus 


(2.2.9)  Q > min  {G['(A+d  JG^'"’(d,  )-G*;(A-d  jG'^'"'(-d  J } , 
■"  0<r<m  n I n I n ^ ^ 


By  Lemma  2.2.1,  we  get 


(2.2.10)  inf  ag(l)  = min{G[;'"’(d^)-6'''"'(-d2),GjJ(A+d^)G'^'"’(d^)  - 

k-m/ 


G';'(A-d2)G"-'"(-d2) 


This  completes  the  proof. 


Theorem  2.2.1 . If  the  family  of  density  functions  {gp(t-e);  Q k h-^} 
has  montone  likelihood  ratio,  for  every  n 1 , then  for  arbitrary 
fixed  d.|,d2  satisfying  d.|+d2  > a, 

(2.2.11)  inf  a,,(l)  = minfG^(d^)-G^(-d2),  Gj^(d^  )-gJ^( -d2) , 

Ci  - 

G[j‘"'(d^  )GjJ(A+d^  )-G^*'"(-d2)G’^(A-d2) , m=l , . . . ,k-l  I . 
k-1 

Proof.  Since  i'  = L ; , by  Lemma  2.2.2,  we  have 
m=0 

/ooio\  * £■  t * r/\k*'mr,k”B  »k“mpmr,k — mp.m  ■, , 

(2.2.12)  inf  -I  (1)  - min  'min[A  -B  , A C -B  DJI 

,/  - i'- 


4 
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where  A = B = G^(-d2),  C = G^(/>+  dj)  and  D = G^^lA-d^). 

By  Lemma  2.2.1,  min  = min[A-B,  A^-b'^].  From 

0<m;_k-l 

(2.2.12) ,  we  get 

(2.2.13)  inf  x,(l)  = min{A-B,A‘^-B^A'^"'"c"’-B‘^'V,  m=l  ,2, . . . ,k- 1 . 
This  completes  the  proof. 


Corol lary  2.2.1.  Let  the  family  of  density  functions  {g^(t-o): 
■ c ’'i  have  a montone  likelihood  ratio,  for  every  n ^ 1.  For 
arbitrary  fixed  dp  d^  satisfying  g(-d.,)  g(d-|)  and  cl^+d^ 

then 

(2.2.14)  inf  <.,(1)  = min{G^(d^)-G^(-d2),G|j(dp-G[;(-d2), 


Proof.  By  Theorem  2.2.1,  we  get 

(2.2.15)  inf  ^t^^(l)  = min{A-B,A‘^-B'^,A‘^‘'"c'^-B'^''^D'^,  m=l,...,k-l} 


where  A = Gn(di)i  5 = G^C-d^),  C = G^(a+o^)  and  D = G^( 
Consider 


(2.2.16)  f(m)  = - b'^''"d'^ 


rtk/C\n)  r,k/D\m  ni  i i 

= A (^)  - B (g)  , m = 0,1 ,. . . ,k-l  . 


Let 


i>(x) 

d'k(x) 

dx 


9n(x+d^) 

G^lx-d^) 


fG  (x-d-y)-G^(x+d, ) 
n 2'  n'  r 


^x+dlT' 
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g(x-d2) 

since  > -d^ , ^ is  nondecreasing  in  x > 

then  ^ 0 for  x > 0.  Thus  oi'-)  1 v(0) 

By  Lemma  2.2.1  and  1 ^ then 


0 and  g(-d2) 

A 

B- 


g(d^ ) , 


(2.2.17)  min  f(m)  = niin{f(0),  f(k-l)}. 

O^m-'k-l 

Hence  from  (2.2.15)  and  (2.2.17), 

inf  ajl)  = min{A-B,  A*^-B^  f(0),  f(k-l)} 

k k k-i  k-1 
= min{A-B,  A -B'^,  AC*^  -BD 

This  completes  the  proof  of  the  corollary. 

Now  we  wish  to  discuss  the  probability  that  contains  all 
good  populations  and  the  infimum  of  this  probability  over 

Lemma  2.2.3.  If  the  family  of  density  functions  'g^(t-  ■): 

■ C ■ 1 ; has  a monotone  likelihood  ratio,  for  every  n 1,  then  for 
arbitrary  fixed  d-j  ,d2  satisfying  d-|+d2  > A, 

(2.2.18)  inf  a (2)  ^ min  min  {[G^(d.|  (-A+d^ ) 

;7  - O-  r'k-m-2  0<i'm 

m 


Proof.  For  o e 


•■m’ 


(2.2.19)  ^^,(2)  = t 


*^^[k]'^l  ’"[k]^'^2^’  ^ ' 


" ^'''[m+1]  - ’^[k]''^r  Ik]  - ^[kf'^2' 
k k 

= ^n("[m+l]-''[jf'^l^  - ^n("[k]-°[j]-^2) 

k k 

- :■  G (•.  -■'  ,1^dJ  - ;i  G„(.s.-dJ 

n'  j m+l  r n'  j 2 
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where 


j " Tk]  - -[j]’  J 


1 , . . . ,k  and 


k k-1  ■ ■ ■ m+1  fii  ■ ■ ■ 1 ■ 

Since  i.  (2]  is  nonincreasinq  in  where  ■-^^2  '■  ] 

I (?)  • Q where 


(2.2.20)  Q - A -•  G f ■ .-r+d,)-B  G ( ■ .-d.), 

j=m+2  " J ' ^ 


j=m+2 


A = n G (■'.■-;:+d-,  )G  (djG  (-A+dJ  and 
j_-]  n'  j r i r 

m 


For  i = m+2 , . . . ,k-1 , 


k-1  k-1  g (^.-dp) 


j2_  = 


g„(  1-^2^ 

Again,  since  d,  > ;.-d,  and  — n — ,---,- -v-  is  nondecreasing  in  2,., 

- 1 '’r  ^ 

i = m+2, . . . ,k-l , then  Q is  greater  than  or  equal  to  the  right  hand 
side  of  equation  (2.2.20)  with  either  2^.  = 0 or  = a,  i = m+2, 
...,  k-1. 

Let  H,  = ^ ^m+’ ’^m+?’ ■ • • ’ ’'k-1  ^ ^ ■’  'k-1  'ni+1  ' ‘ 


Hence, 

(2.2.21)  inf  a (2)  - min  CAG^(-/+d,  )G|^''""^'^(d^  )-8G|;^(-d2) 
H,  - 07r<k-m-2 


Let  Ho  - 


. 2,  \ 

"m  ■ ■ ■ ■'  1 ’ ' 


In  a manner  similar  as  above,  for  r - 0, 1 , . . . ,k-m-2 , we  can  show 


that 
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(2.2.22) 


inf[AG[J(-A+d^)G[^""'"^"'"{d^)-BG[;(-d2)G 

^2 

= nn-n  f6;;"^-A.d,)[G 
0<;.^iii  ^ I n I 


k-m-2-r/ . , 
n 

ir  I.  A Mk-ni-l-r+£ 

-LG^lA-d^)  J 


From  (2.2.21)  and  (2.2.22),  wo  have  our  lemma. 
k-1 

Since  sf  U s;  , by  Lemma  2.2.3,  we  can  have  the  following  theorem. 

0 1 II 


Theorem  2.2.2.  If  the  family  of  density  functions  {gp(t-o):  o C;  A.|) 
has  a monotoru;  likelihood  ratio,  for  every  n ^ 1,  tlien  for  arbitrary 
fixed  dpd2  satisfying  d^+d2  > A, 

(2.2.23J  inf  '<,  (2)  = min  min  min  l[G  (d,)]'^^*'  ^ 

\i  - 0-^’m-:k-l  0--r<k-m-2  0<y,<iii 


Lemma  2.2.4.  Suppose  that  (i)  9p(t)  = 9p('t)  > 0 a11  t and 

(ii)  the  family  of  density  functions  {g^(t-a):  0 L A^l  has  a 
monotone  likelihood  ratio,  for  every  n jy  1 . For  arbitrary  fixed 
d-j,d2  satisfying  d-j  = ^2  = d > A, 

(2. 2. 2^1)  inf  («|,(2)  = mi B-c'''"'"' ^ U,  a'^’^B-c'^’^ D,  Ab'^""’’^- 
^■'■m  - 

^pk-m-1  ^m+lpk-m-1  ^m+lj^k-m-lj 

where  A = G|^(d-|),  B = G|^(-A+d|),  C = 6^(/'.-d2)  and  D = Gj.(-d2). 


Proof.  By  Lemma  2.2.3,  we  have 


z',  o or\  ■ r • • , a k+ I’.-iii- r - 1 „r +1  ^ k+ p.-iii- r- 1 p. 1 

(2.2.25)  inf  n„l2)  mm  mm  (A  B -C  D 

O' r-d;  •m-2  O'  v.-.m 


where  A,  B,  C and  D art?  defined  as  above. 
Consider 
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(2.2.26)  Q = min 
Now  we  want  to  show 


0-  r-  k-m-2 
B D 


A L' 


Let  f(x)  - G (x-d)-C/.(.v-d)  for  G x • d. 

f t II 


df(x) 


' gr,(x-d)[l-2G^,(x-d)  ] C,  0 • X ■ d. 


Since  0 • d,  f(')  ■ fi'O)  .jnd  G^it)  ' l-G^(-t),  we  can  see 


that 


B D 


A ■ C' 

By  Lemma  2.2.1  and  from  (2.2.26), 


(2.2.27)  Q :=  min<A^^'""-^B-C^^'‘'”-^D,  A' b"'""' ^ ^ . 
^roni  (2.2.25)  and  (2.2.27), 

(2. 2. 28)  inf  ,,,(2)  --  min'  min 


B-C 


0-  . -m 


'D], 


nnn  LA’ ^ ' ] 


B / nr-i[::a  2.2.1  a n d ( 2 . 2 . 28 ) , we  qe  t 


(2,2,21j  mf  . (2)  --  min  rni nLA*'''"' ' D, A^' ' ^ B-c'^' ^ D] , 


.;„i  1 /,,.k-m-l  .m+1  „k-ni-l  ^m+lr.k-m-1 

iinLAD  -oD  ,A  E -C  D 


This  completes  tiie  proof. 

k-  1 


>1  nee 


, bv  I omma  2.2.4  and  Lemma  2.2.1,  we  can 


obtain  the  follow  inn  result. 


Corollary  2.2.2.  Suppose  that  (i)  (||,(t)  - ^ ^ 


'o.:  'ii)  the  faiiiily  of  (lensity  functions  ‘ q^.|(  t.-- ) : “ c has  a 


I 

J 
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monotone  likelihood  ratio,  for  every  n > 1.  For  arbitrary  fixed 
dpd^  satisfying  d^  = d2  = d > a, 

(2.2.30)  inf  a (2)  - 

fi 

where  A,  B,  C and  0 are  defined  as  in  Lemma  2.2.4. 

Now  we  wisli  to  discuss  the  probability  that  I-j  excludes  all 
the  bad  populations  and  the  infimum  of  this  probability  over  S2. 


Lemma  2.2.5.  For  o G , 

- 111 

(2.2.31)  «q(3)  > l-G^-'"(d^)G(-A+d^) 


Proof.  For  (I  C-  Q , 
m 


(2.2.32)  a^^(3)  = P[0(-.j  Mr  i = l,...,m] 

- - ^[k]'^'^2^^''^°[m]  - "^[k]'^!^ 

k k 

Let  6.  = i = 1,2,. ..,k. 

k k 

(2.2.33)  ni^.j(3)  = ^n^'^j''’'l~'^2^~  ^n^ "'^m^'^1  ^ 


(<a,3(3) 


since  - 


^ 0 for  6.|  > ^2,  by  letting  5.|  = ™, 


then 


(2.2.34)  ^3)  . I - ,«  G„(YV^1> 

0 ^ 


Since  the  right  hand  side  of  (2.3.34)  is  nonincrcasing  in  6.  for 

J 

j = 2,3,...,m-l,  m+l,...,k-l,  and  is  nondecrcasing  in  it  follows 


° '^k-1  Vl  ^ ^ *^2^  ■'^  achieved  at  =...= 

‘^k-1  " '^m  " ^ ^2  "•■•"  Vl  " “■ 


(2.2.35)  a^(3)  > 1-G|;-'"(d^  )G^(d^-A) 

The  proof  is  complete. 

k-1 

Since  n = u by  Lemma  2.2.5,  we  can  easily  lead  to  the 
m=0 

following  theorem. 

Theorem  2.2.3.  For  arbitrary  fixed  d^d2  satisfying  d.j+d2  > 0, 

(2.2.36)  a^O)  I l-G^(d^)G^(-A+d.|). 

Remark  2.2.1.  (i)  The  right  hand  side  of  (2.2.36)  is  independent 

of  d2*  That  is,  for  a given  d.|  and  for  any  fixed  d2  satisfying 

d^+d2  > 0,  Oq(3)  is  bounded  below  by  1 -G^(d.|  )G^(-A+d^ ) . If  the 

confidence  interval  is  defined  as  I.|  = (Tj-j^j-d-,  ,co)  where  d.j  is 

any  given  number,  then  inf  «^(3)  can  be  obtained  as  follows. 

n - 

(2.2.37)  inf  a^(3)  = l-6^(d^ )G^(-A+d^ ) . 

(ii)  If  we  assume  ag(4)  is  the  probability  that  I.|  contains 
at  least  one  good  population  and  also  I.|  excludes  all  bad  popula- 
tions, then  f^g(4)  ,^0^(1)  + «g(3)  - 1.  From  (2.2.11)  and  (2.2.36), 
we  get 

a^(4)  imin{G^(d^)-G^(-d2),G|^(d^)-Gj(-d2).G|;"'"(d^)Gj;(A+d^)  - 
Gjj‘'"(-d2)G'"(A-d2)  ,tn=l , . . . ,k-l  }-G(d^  )G^(-A+d^ ) . 


Similarly,  we  can  find  the  lower  bound  of  =<^(5)  which  is  defined 
as  the  probability  that  I.|  contains  all  good  populations  and 
also  I-j  excludes  all  bad  populations. 

Let  S'  be  the  number  of  non-best  (bad)  populations  that  enter 
Wo  are  interested  in  E(S'),  the  expected  size  (over  fJ)  of 
non-best  populations  that  enter  I-j. 

Theorem  2.2.4. 

(2.2.38)  sup  E(S')  £ (k-l)G  (d^) 

Proof.  For  0 t 

m 

(2.2.39)  E(S)  = P[0...  t 1} 

i = 1 ^ 

m k k 

k k 

Gn(-d2)G[j 

im{G^(d^)G|;-'’’-^(A+di)G';(o^k]-°[lfdl^-Gn(-d2)- 

'^n  ^[^°[l]"®[k]'^2^^ 

< 111  G^(d^)G^;"’-''(A+d,). 
k-1 

since  = u then 
m=0 

(2. 2. TO)  sup  E(S')  - maxim  G^(d^  )G^^'"''^  ( A'l-d^ ) iii=0,l , . . . ,k-l  1 
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In  the  following,  we  deal  with  special  values  of  d-j  ,d2  and  a. 

Case  ( i ) : Arbitrary  but  fixed  d^  ,d^  satisfying  d-j+d^  - 0 and 

A = 0.  Since  A = 0,  the  only  good  population  is  the  one  associated 
with  Applying  the  earlier  results, 

(2.2.41)  inf  xjl)  --  in"  , (2)  - siin  , Gjd^  )-G^(-c^)  ,G|;(d^  )-G^( -d^) , 
a result  obtained  by  DudewicA  and  Tong  [26]-  Also,  in  this  case 
{2-.2A1)  .{3)  1-G^"(d,). 

Case  (ii):  One  side  (lower)  confidence  interval.  Arbitrary  but 

fixed  d^  and  d...  - ' with  ■ ■ 0.  Our  confidence  interval  is  now 


(2.2,43) 

i n f ■ ( , 

(1)  ^ 

(2.2.44) 

inf  -t 

(2)  - 

wV' 

-■■■"dl 

) and 

(2.2.45) 

inf  1 

(3)  - 

l-'^  ^d,)G  ( 

-'"dl 

). 

Case  (iii) 

: One 

side 

(upper)  con 

f iden 

ce  interval.  Arbi trary 

fixed  d2, 

dl  ^ 

and 

■ • 0.  Now 

our  confidence  interval  is 

defined  as 

M = 

(-■, 

T[^]+d,),  wo 

can 

snow  that 

(2.2.45) 

inf  : 

= l-Ci-d.), 

(2.2.46) 

inf  t 

,(?) 

- i-i.i  (.  -d 

n 

,)G(- 

ri„;  and 

(2.2.47) 

inf  a 

,.(3) 

- 0. 

ii 
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’.3.  Result  for  the  scrile  re »''ii"e'''er  cere. 


In  this  section,  we  uscuiiif-  j ■,  a scale  parame^’er  for  T ■ , 

: " l,...,k.  It  is  assu'sed  that  there  i s no  a priori  knowledge 
about  tne 

.G  ('  1 for  t 0 

I2.3J)  G„(t,  ,)  .{  < 

0 for  t ■ 0, 

for  ■ 0,  k = I , . . . ,k. 

Tor  arbitrary  fixed  a,  n --.a tr. lying  a 1 and  b 1,  if  the 
confidence  interval  is  define:!  a.. 


(2.3.2) 


^2  " 'a  " ^[k]^ 


wnere  Tr,  n = irax  T • . 

lj;k  ^ 

Let  be  a given  positive  (onstarit,  -1.  We  say  population  is 


' •?  ■?  "5  1 

'y  • > . a ; 


good  i f 


bad  if 


1-,  k ' 

For  given  P and  I,  lef  (1)  denote  tfic  probability 

that  It  contains  at  least  one  good  pojnilation.  let  . (2)  denote  tlv' 

probability  that  [.,  contain,  ail  OMod  population  and  let  : (3)  denote 
the  probability  that  I exclii.aes  all  the  t;ad  populations.  Using  the 
saiae  argunents  as  in  Sort  ion  2 . ? , v»e  outain  the  following  results. 
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Theorem  ?.  3 . 1 . If  the  family  of  density  functions  (g^(-;):  o (.‘.2'! 
has  monotone  likelihood  rotio,  for  every  n 1 , then  for  arbitrary 

fixed  a,  b satisfying  a ' 1,  b 1 and  ab  ■ , then 

(2.3.4)  inf  '.  (1)  ^ „in  G^(a)  - 

(a)u  (a/,j-G^  (g^)G^(-g-) , m=l , . . . ,k-l ; . 

Coro! lary  2.3.1.  Let  the  family  of  density  functions 

■■  V has  monotone  likelihood  r(3tu),  for  every  n > 1.  For 

arbitrary  fixed  a,b  satisf/inq  a • 1,  b ’1,0  (r)  > g (a)  and  ab  > A 

(2  3.5)  inf  , (1)  ^ min  G^(a}-G,,(j^),Gj;(a)-G|;(i),  G^(a  )g|;' ' (aA  )- 

'“.:"i1y  of  density  functions  {g|^(“):  '■  t • 
has  a 'lionotono  likelihood  ratio,  for  every  ■'  • 1,  then  for  arbitrary 
fixed  a,b  satisfyirig  a 1,  i;  ■ ■ and  ai> 


(i^.3.6)  inf  , (2)  = i;:in  min  aiin  'fG  (a)] 

- n 


k+ . -m- r- 1 


f -1  O r k-m-? 


.-r ’■  I ,a 


Theorem  2.3.3.  For  arbitrary  fixed  a i and  b ■ 1, 


(2.3.7) 


(3)  ■ l-G„(a)G  {^) 


Lets'  be  the  number  of  non-best  (bad)populations  that  enter 
Also  we  can  obtain  the  followmq  result. 


T h e 0 rem  2 . 3.4. 

(2.3.8)  sup  E(S')  ■ (k-l)G^(a). 

In  the  following,  we  deal  with  some  special  values  of  d.| , 
d^  and  a. 

Case  (i ) : Arbitrary  fixed  a,b  satisfying  a 1 and  b ^ 1 and  = 1 . 

Since  \ - 1 the  only  good  population  is  the  one  associated  with 
Applying  the  earlier  results, 

(2.3.9)  inf  ;,..(1)  = inf  - ,,(2)  '"i  n- G^( a ) -G^ (-^)  ,G^ (a ) -G^ (g-)  ; , 

a result  obtained  by  Saxcna  and  Tong  [66].  Also,  in  this  case 

(2.3.10)  :•  ,(3)  ^ 1-G^Ia). 

Case  (ii):  One  side  (lower)  confidence  interval.  Arbitrary 

but  fixed  a • 1 and  b = > with  ■ 1.  Our  confidence  interval 
is  defined  as  lo  ‘ (a  Trp]’'i‘  that 

(2.3.11)  inf  , A\)  - g"^(i), 

(2.3.12)  inf  ;-,,(2)  G^,(a  )gJ^' ' ([' ) and 

(2.3.13)  inf  f,,(3)  - l-Gj^(a)G^^n  . 

Caje^X-iij):  One  side  (upper)  confidence  interval.  For  arbitrary 

fixed  b a = ” and  , 1.  New  ou^  n:H\fu!encc  interval  is  defined  as 
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We  can  show  that 


(2.3.14) 


inf  : (1)  - 


1-G  ({-), 

nij' 


(2.3.15) 


inf  ; (2)  - 


1-G  and 

n ' b n b 


(2.3.16) 


inf  , . (3J 


2.4.  Some  examples. 


We  shall  discuss  the  above  problem  for  the  means  of  normal 

populations  and  for  the  scale  parameter  of  the  gamma  population. 

Let  "i  , be  k normal  populations  with  unknown  means 

2 

■'■■j , V2,  • ■ . . and  a co;imon  known  variance  . Let  X.j  be  the  sample 
mean  from  based  on  a sample  of  size  n.  Let 


(2.4.1)  T . = X . and  T 


1 1 


[k]  - 


max  7.. 


Let  G|^(x,n^.)  be  the  c.d.f.  of  from  i = l,...,k.  In  this 


case,  G^(x,-.)  = G^(x-'. ) and 


(2.4.2) 


G (x)  ^ :(-  x) 
n 


where  ^ is  the  standard  normal  distribution  function.  The 
confidence  interval  I|  is  of  the  form  Tj-|^j  + d„), 

where  d.j+d^  > . !n  this  case,  applyinn  Theorem  2.2.1,  we  get 

(2.4.3)  inf  , (1)  - nhn-A-B,  a'^-b'",  , nc  1 , 2 k-1 


where 


-.n  d. 


. n(  +d. 


A - :( f),  B = M.  C = and  D - : ( 


1 :i(  -d.,1 
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For  a given  P*,  A > 0,  > 0,  > 0 and  d^+d^  > A,  since 

1 im  A - 1,  lini  B “ 0,  1 im  C - 1 and  1 iin  D = 0,  there  exists  a 
n- n n-^-"  n-'-“ 

smallest  integer  N such  that  for  every  n ^ N,  we  have  inf  a (1)  > P*. 

n - “ 

Let  X|  be  the  unique  solution  of  '»(d^x)-'?(-d2x)  = P*.  Let  X2  be 

k k 

the  unique  solution  of  ■!'  (d-jx)-':'  (-d2x)  ^ P*.  Let  be  the  unique 
sol ution  of 

*^’'""^(d^x)t'^((A+d^)x)-'!>'^~'"(-d2x)<r."’((A-d2)x)=  P*,  m=l  ,2, . . . ,k-l 
Let 

N.  = [n"x^],  i = 1,2  and  M.  [a^y^],  i = l,...,k-l, 

where  [x]  denotes  the  smallest  integer  j x.  It  follows  that 

the  smallest  sample  size  N required  to  satisfy  inf  « (1)  > P*  is 

- 

given  by 

(2.4.4)  N = iiiiniN.,  M.,  i=l,2,  j = l , . . . ,k-l } . 

' sJ 

Consider  the  probability  that  I-|  contains  all  good  populations. 

For  given  P*,  a > 0,  d.|  :•  a and  d2  '•  0,  then 

lim  Gn(d-|)  = 1,  lim  G^(-A+d.j)  = 1 and  lim  G^(-d2)  = 0. 

By  Theorem  2,2.2,  we  can  always  find  the  smallest  sample  size  N 
which  satisfies 

(2.4.5)  inf  r.^^(2)  > P*. 

Consider  the  probability  that  I.|  excludes  all  good  populations, 
for  given  P*,  A > 0,  d^ , dp  such  that  d^+dp  > 0 and  d^  A,  then 

1 im  G(-A+di ) = 0. 
n 
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By  Theorem  2.2.3,  we  can  always  find  the  smallest  sample  size  N 
required  to  satisfy 

(2.4.6)  inf  d*  . 

Now  let  us  considei"  tnr  uroblei?:  in  relation  to  the  scale  ; 

i 

parameters  of  tne  population;.,  uot  be  k gamma 

distributions.  The  distribution  associated  with  has  the 
density  function  g(x,-,,  .)  where 


(2.4.7) 

g ( X , It , a . ) 


for  X > 0 


otherwise. 


We  assume  that  the  shape  parameter  ^ is  the  same  for  all  populations 
and  -t  is  known.  Let  T.j  be  the  sample  sum  from  based  on  n observa- 
tions and  let  Tr,-j  - max  T-.  Then  the  confidence  interval  is 

1-i-k  ’ 


given  by 


(2.4.8;  b .-.b 

Let  G^(x,"^. ) be  the  c.d.f,  of  T.  from  i “ l,..,,k. 
In  this  case,  G^(x,-.i.j)  ^ 


(2.4.9)  distribution 

with  scale  parameter  1 and  '.iiape  p-iraibeter  n-,..  Since  the  condition 

in  Theorem  2.3.1  holds,  we  get 


12.4. in)  inf  (1) 


min ; n( a ) -G ; , 


(a;- 


G (a)G  (a',)-G  (|JG  (|^ ) ,m=l  , . . . , k- 1 t . 


For  given  P*,  ; 1,  a 


JWUUIIWUIILI 


1 im  G [r]  = 0 and  1 ur  G (a  ) 
n b n ' 

n n - 

N such  that  for  every  n N, 


I with  ab  • , since  inn  G 'ai-', 

n 

n -■ 

= 1,  there  exists  a smallest  integer 
we  have  i nf  / ( 1 ) > P* . 


Let  be  the  smallest  positive  integer  such  that  Gp(3)- 

G^(^)  P*.  Let  M2  be  the  smallest  positive  integer  such  that 

^ k 1 

G^^(a)  - G^^(^-;  ■ P*.  Let  M,n  be  the  smallest  positive  integer  such 

tnat  G|^''’'(a)  G'''(a.')  - ) , m - 1 ,k-l . . 

It  follows  that  the  smallest  saisple  size  N required  to  satisfy 
inf  I (1)  P*  is  given  by 


(2.4.11)  M = min'N.,  i = l,2,  j = l ,2, . . . ,k-l  } . 

Consider  the  probabiiit.  that  I-,  contains  all  good  populations. 
For  a given  P*,  ' • 1,  a • ' and  b ■ 1,  then 

lim  G (a)  - 1,  lim  G i ■'  1 and  lim  (fg)  = 0. 

n >■  ' n - ' n - ■ “ 

By  Theorem  2.3.2,  vie  can  always  find  tne  smallest  sample  size  M 

required  to  satisfy  inf  (';)  - b* 

Consider  the  nrobabilit/  that  I,  excludes  all  good  populations. 

Por  given  P*,  1,  b ' •tsi  .1  1 .vi  th  a , then  lim  G (^)  = 0. 

n ' ‘ 

fly  Tneorem  2.3.3,  we  can.  .-il.vays  riou  th(>  simillest  sample  size  .M 
required  to  satisfy  inf  :•  ( 3 1 P*. 


Remark  2.3.1.  (i)  IP  the  sample  si.’p  M is  preassiqned  and  P*  and 
■ arf?  given,  then  the  confidetue  interval  ran  always  he  found 
fron  (2.2.11),  (2.2..P3)  ,ind  .IG)  with  Gj.|V.r)  as  defined  in 

^2.4.2)  so  that  inf  g , ( i ) ■ P * , i 1,2,3. 
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Similarly,  the  confidence  interval  can  always  be  obtained 
from  (?.3.4),  (2.3.6)  and  (2.3.7)  with  0|,(x)  as  defined  in 
(2.4.9)  so  that  inf  . .. (i)  ■ P*,  i = 1,2,3. 

(ii)  In  location  parameter  case,  if  tf:»  population  with  the 
smallest  ■!- value  is  c<tlled  the  best  population.  Let  '.  be  given, 
we  say  population  is 

(2.4.12)  good  if  • . ' t ] - . 

bad  if  . - * 

L ' J 

where  t-rii  ‘ 

l;'J:'_k 

Then  the  confidence  interval  is  given  sy 


(2.4,13) 


1 


(T 


[1]' 


where  Tr 


[1] 


min  T.,.  . 
1 i ' k 


The  calculation  of  the  coverage  prooabi  1 i t i t-s  g,(i),  i - 1,2,3, 
and  their  infima  can  be  nandled  sinni’arly  as  in  Section  2.2. 
The  scale  parameter  case  c n i>e  nandlei  in  a similar  manner. 
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CHAPTER  III 

SELECTION  PROCEDURES  IN  TERMS  OF 
MAJORIZATION  AND  WEAK  MAJORIZATION 

3.1.  Introduction 

Over  the  last  fifty  years,  majorization  and  Schur  functions 
have  been  applied  to  develop  many  useful  inequalities  in  many 
branches  of  mathematics.  Some  of  the  references  are:  Schur  [67], 

Ostrowski  [57  ],  Hardy,  Littlewood  and  Polya  [45],  Beckenbach  and 
Bellman  [17]»  Marshall  and  Proschan  [52]  and  Marshall  and 
01  kin  [50  ].  In  recent  times  the  techniques  of  majorization  and 

Schur  functions  have  been  applied  to  probability  and  statistics 
(see,  for  example,  Marshall,  Olkin  and  Proschan  [51],  Rinott  [62], 
Proschan  and  Sethuraman  [61],  Nevius,  Proschan  and  Sethuraman 
[55]  and  Gupta  and  Wong  [44]). 

Recently,  Nevius,  Proschan  and  Sethuraman  [56]  introduced 
a stochastic  version  of  weak  majorization.  They  have  discussed 
some  properties  and  made  some  applications. 

Gupta  [32]  defined  a class  of  selection  procedures  and 
considered  some  of  its  properties.  Some  additional  results 
concerning  the  properties  of  this  class  of  procedures  were 
obtained  by  Gupta  and  Panctiapakesart  [36],  Gupta  and 
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Panchapkesan  [37]  defined  a class  of  selection  procedures  which 
is  a natural  generalization  of  the  class  considered  by  Gupta  [32]. 
In  that  paper,  they  obtained  a sufficient  condition  for  the 
probabi!  i ty  of  a correct  'ioloction  to  be  nondecreasing  in  a when 


A,  where  a^-  is  the  parameter  associated  with  ir^. , 
They  also  obtained  the  supremum  of  the  expected 


size  of  the  selected  subset  and  showed  that  if  the  sufficient 


condition  holds,  it  takes  place  when  the  A^-'s  are  equal. 

In  this  chapter,  we  order  the  parameter  space  by  means  of 


majorization  or  weak  majorization  and  propose  some  selection 


procedures  when  the  parameter  a^  associated  with  , i=l,...,k 
is  a vector.  Section  3.2  defines  a class  of  procedures  for 
selecting  the  population  associated  with  vector  a^j^j  (see 
Section  3.2  for  definitions).  A sufficient  condition  is  obtained 
for  the  infimum  of  the  probabi  1 ity 'of  a cor»'ect  selection  to  be 


Schur-convex  in  A.  Also  another  sufficient  condition  is  obtained 


for  the  supremum  of  the  expected  size  of  the  selected  subset  to 
take  place  when  a.|  =...=  Aj^.  Some  special  cases  of  interest  are 
discussed.  In  Section  3.3,  a sufficient  condition  is  obtained 
for  the  same  infimum  of  the  probability  of  a correct  selection 
to  be  nondecreasing  and  Schur-convex  in  A.  Section  3.4  defines 
a class  of  procedure  for  the  selection  of  the  population  with 
vector  Some  properties  of  the  selection  procedure  are 

briefly  discussed.  Section  3.5  and  3.6  deal  with  selection  procedures 
for  multivariate  normal  distributions  in  terms  of  majorization 
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and  weak  iiiajorization  of  mean  vectors.  Various  cases  corresponding 
to  the  known  or  unknown  common  covariance  matrix  Z are  studied. 
Properties  of  those  selection  procedures  are  also  established. 

3.2  A class  of  selection  procedures  R^^(and  some  properties) 
for  vector-valued  in  terms  of  majorization 

First  we  give  the  definitions  of  majorization,  Schur-concave 
and  Schur-convex  functions. 

Definition  3.2.1 . A vector  a = (ap...,ap)  is  said  to  majorize 

a vector  b - (b^ , . . . ,bp) , if  a^  j;.-  • -il  t>-|  bp,  and 

r r p p 

y a.  > y b.,  r = l,...,p-l,  while  )’  a.  = ) b.;  we  write, 

i=l  ^ ■ i=l  ^ i=l  ^ i=l  ^ 

a > b. 
m 

The  above  definition  is  according  to  Beckenbach  and  Bellman 
[ 17]  and  differs  slightly  from  that  of  Hardy,  Littlewood  and 
Polya  [45]. 

Definition  3.2.2.  H is  called  a Schur-concave  (Schur-convex) 

function  if  H(o)  H(b)  (H{a)  2,h(b))  whenever  a > b.  A Schur 

~ m ~ 

function  is  a function  which  in  either  Scfiur-concave  or  Schur- 
convex  . 

From  above  definitions,  we  know  that  majorization  is  a partial 
ordering  in  the  p-dimensional  tuclidean  space  and  Sciiur  functions 
are  functions  that  are  monotone  wi  tti  respect  to  this  partial 
ordering.  Mow  we  state  the  following  theorem. 
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Theorem  3. 2. 1 . (Ostrowski  [57]) 

Assume  H is  defined  for  2^  and  first  partial 

derivatives  of  H exist.  Then  H(Z)  - H(Z')  for  all  Z Z'  if 

m 

and  only  if 

(3.2.1  ) ^ - y-  '■  0,  for  i ■-  j , i , j = l ,. . . ,p. 

i j 

Let  - be  k populations.  Let  A be  an  interval  on 

the  real  line.  Associated  with  ■! . is  a vector  X.  = (X.t,...,X.  ), 
i 1 , . . . ,k  where  X . • has  density  ;(■•■, x),  i = l,...,k,  j = l,...,p. 

I J 1 J 

We  assume  that  C and  '.i  •...  • ' -p,  i = l,...,k.  We  say  that 

1 j -1  ' ir  ’ ip'  -j  ^ ji  ’ ’ jp^ 

It  is  assumed  that  among  the  k given  populations,  there  always 
exists  one  population  which  is  better  than  others.  This  is  equiva- 
lent to  saying  that  there  exists  a such  that  for 

all  i = l,...,k.  This  population  is  called  the  best  population. 

If  there  are  more  than  one  "best"  population,  then  we  assume  that 
one  of  them  is  tagged  as  the  best.  Based  on  one  observation 
^il’.’.’Xip  population  we  construct  a suitable  statistic 

T.  - g(X.p...,X.p),  i = 1 k.  Let  F,  (x)  = P[T-  < x]  be  the 

distribution  function  of  T-  from  population  Based  on  the 

values  of  T^-  from  i - l,...,k,  we  wish  to  define  a class  of 

procedures  for  selecting  a non-empty  subset  of  the  k populations 
such  that  the  probability  is  at  least  P*  (p  < P*  '■  i ) that  the 
population  associated  with 


is  included  in  the  selected  subset. 


Let 

(3.2.2)  E = (A  = (Ap...,Ap):  A^  C A.  i = ',..., p} 

and 

(3.2.3)  Q = (w  = (A^,...,Aj^):  A.J  tE,  1 = 1,. ..,k  and  there  ex1  sts 

some  1 such  that  A.  > A.  for  each  j}. 

' m 

We  wish  to  define  a selection  rule  R such  that 

(3.2.4)  inf  P[CSjR]  > P*. 

Let  h * h c <;  [1 ,~) , d € [0,®)  be  a function  defined  on 
the  real  line  satisfying  the  following  properties:  for  every  real 

X, 

(3.2.5)  (i)  h^^j(x)  1 X 

(ii ) h^  ^q(x)  = X 

(iii)  h .(x)  is  continuous  in  c and  d 

C 9 Q 

(iv)  h^  ,(x)  t 00  as  d ->  «>  and/or  xh  .(x)  t <»  as  c xfO. 

C ^ Q C > U 

Some  functions  satisfying  these  properties  that  will  be  of  interest 
are  cx,  x+d  and  cx+d.  Now,  we  define  a class  of  procedures 
R^  as  follows. 

R^:  Select  population  if.  if  and  only  if 

(3.2.6)  h(T.)  > max  T . 

l£r<k 

Because  of  (3.2.5)- (i),  the  procedure  will  always  select  a 
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non-empty  subset.  Let  be  the  random  variable  associated 
with  'r.-i  and  let  F,  be  the  cumulative  distribution  function 

-LiJ 

(c.d.f.)  of  i = l,...,k,  we  have 

(3.2.7)  P[CSjR^]  = P[h(T(^.)  i^l,...,k-l] 

- k-1 

= / r,  (h(x))dF,  (x). 

-[i]  -[1] 

We  now  assume  that  the  distributions  are  Schur-conca ve  in  x, 

i.e.,  for  X > X',  F,  and  F.,  are  distinct  and 
m ' - 

(3.2.8)  F,(x)  F^  ,(x),  for  all  - . 

Then 

P[CS.RJ  > / F*"'^  (n(x))dF  (x). 

^ -[k]  -[k] 

Hence 

(3.2.9)  inf  P[C$iR^]  - inf  .(_;  c,d,k) 

5.'  ' 

where 

(3.2.10)  v'(7;c,d,k)  - / F^'*^h(x)  )dF  Jx) , X tE. 

In  the  fol lowing  theorem,  we  are  interested  in  a sufficient 

condition  for  the  Schur-convexi ty  of  E,(T,a)  where  a(x,x)  is  some 

real-valued  function,  x = (x,,...,\  ) F and  T=g(X.  ,...,X.  ) is  a func- 

P 'l  'p 

tion  of  X ,...,X  . Let  F (x)  = P[T  ■ <]. 

At  .\  \ 

1 p 

Theorem  3.2.2.  Let  'Fjx),  '■  C E;  be  a family  of  continuous 
distributions  on  the  real  line  such  tnat  F (x)  is  a differentiable 
function  in  x and  ■.  and  let  :(x,'-)  be  a bounded  real-valued  and 
differentiable  function  in  x and  • , i l,...,p.  Then  E,(T,0 
is  Schur-convex  in  . c E provided  that 


r 
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(3.2.11) 


F,(x) 


(~--^-)Fjx)  ( '-).(x,^) 

I '1  j - i j 


0,  for  i-j  ,i , 


ProQ-f.  We  argue  alonq  the  lines  of  the  proof  of  Theorem  2.1  of 
Gupta  and  Panchapkesan  [37].  We  assume  that  F (x),  >.  e F has  the 

A ~ 

support  I.  Let 

(3.2.12)  A(a)  = L,(T,>)  = / ,(x,\)dF.(x). 

1 ■ - 

Consider  u-,,  u,  CE,  where  u • - (u,- 1 , . . . ,u  • ) , i = 1,2, 


^1  - ^2 


(3.2.13) 

(3.2.14) 
and 

(3.2.15) 


ip  U2  1-  , «M..I  t-  U^.  - 1 ’ ■ • • >-.jpi 

11-jp’  i - 1 .2  and  assume  u^  = u^.  Define 


^2^yi ’y2^ " ’yi ^^^u. * 


3(ypu->)  ~ "I" A 2 ( u 1 , y 2 ^ ’ 

By  integrating  A^(upU2)  Gupta  and  Panchapkesan  [37]), 

(3.2.16)  PAupU,)  ■-  a tern'  indonendert  of  Uj 

+ /[.(x,u,)f  (x)- . ' (x,u,)F  (x)]dx 

T 'Jo  - . . U T 


I ' y?  ■'  -1 


f^or  i ' j,  we  have 
(3.2.17)  ( - 


i "'1 


F (x)  - F 

U dX  U 

(x). 

) ('  t (x)[ 

■v(x,u^  ) 

-,'(x  ,U^ 

i 

'"ll 

'F  (x) 

■F  (x) 

yi 

- .‘(x,Lj2. )[ 

-.1.  ] 

;Up 

W 
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Wc  can  nhow  that,  for  . 


(3.2.18)  E(,\,0  = 


- - 

■'  -■  2i=yr- 


where  A = (,\.j , . . . ,Ap)  and  ■ 
From  (3.2. 18)  we  get 


(3.2.19)  B(a,a)  = 2-  E(u,,u,; 

- - -u,.  -i  -2  u,-u. 


and 

(3.2.20) 


'li 


-1  -2 


TT^  B(a.a)  = 2 Bfa,  ,u,)  , . 

>A.  .u,.  '-1 ’-2'  u,  = u.,=A 

J » J - I - 1.  “ 


From  (3.2.19)  and  (3.2.20),  we  get 


(3.2.21)  (^  - 4-)B(^,'•)  = 2(^-  - -t_)B(upuJi  ... 

1A.  ,u^.  ,U^j  -1  -2  'UyU^-X 

If 


(3.2.22) 


3X  U2 


(-f- 

aui-  ju^j 


F,  (x) 


AX 


> 0,  Vi  < j, 

i ,j  = 1 , . .p, 


then  from  (3.2.17), 


(3.2.23)  — )B(u,  ,Up)  ' 0 for  i • j , i ,j-l  , . . . ,p. 

Hence  if  (3.2.22)  holds,  we  have  from  (3.2.21) 


(3.2.24) 
Note  that 


■7-'- r'-)B(',A)  ■ 0 for  i < j,  i,j  = l,...,p. 

'1  'j 


(3.2.25) 


B(a,'.)  --  2 A(A). 


Applying  Theorem  3.2.1  and  from  (3.2.24), 
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it  follows  that  B(a,a)  in  Schur-convex  in  \ if  (3.2.22) 

holds.  Since  = ^2'  Schur-convex  in  A if  (3.2.11) 

holds. 

This  completes  the  proof  of  the  theorem. 

In  some  cases,  we  will  be  dealing  with  the  function 
such  that  it  satisfies 

(3.2.26)  ^ v(x,A)  v(x,A), 

- ''-'p 

then  {j—-  = 0,  i ■;  j , i ,j  = 1 ,. . . ,p.  Hence 

■’  'i  j 

we  have  the  following  result. 

Corollary  3.2.1 . If  satisfies  (3.2.26),  then  E'C;(T,a)  is 

Schur-convex  in  A if 

(3.2.27)  a(x,a)  ^0,  i < j,  i,  j-l,...,p. 

’ " j - 

Corollary  3.2.2.  If  F^(x)  is  Schur-concave  in  ACE,  v(x,a) 
satisfies  (3.2.26)  and  is  nondecreasing  in  x,  then  Eii(T,,\) 

is  Schur-convex  in  A C E. 

Proof.  Since  F^(x)  is  Schur-concave  in  ^ C E, 

(~  - “)F^(x)  <0,  for  i ■ j , i , j-1  , . . . ,p. 

Also  ii'(x,A)  is  nondecreasing  in  v,  then  -J-  i,  (x,0  - 0. 

Hence  by  Corollary  3.2.1,  the  result  follows. 

As  a special  case,  if  :i  (x,  ) (x),  then  we  have  the  following 

result . 


i 
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■O.QIPll^T'y  3.2.3.  If  F^(x)  is  Schur-concave  in  x and  t|j(x)  is 
nondecreasing  in  x,  then  Ei^(T)  is  Schur-convex  in  x. 


Theorem  3.2.3.  For  the  procedure  let  ’'F^(x),  A £ E}  be  a 

family  of  continuous  distributions  on  the  real  line  such  that 
F,(x)  and  F^(h(x))  are  differentiable  I'unctions  in  x and  x., 
i ~ Then  ijj(x;c,d,k)  as  defined  in  (3.2.10)  is  Schur- 

convex  in  X G E provided  that 


(3.2.28)  f^(x)(-^  - ^)F,(h(x))-h-(x)f  (h(x))(4- 

7 • j - ■"S' 

for  i < j,  i,j  = 1, 

where  f^(x)  = ■—  F^(x)  and  h'(x)  = h(x). 


J - 
• • >P, 


0, 


Proof.  By  letting  ij(x,>)  = F^'^(h(x))  and  using  Theorem  3.2.2. 
yields  the  proof. 

We  note  that 


(3.2.29) 


(3.2.30) 


=•••=  ij-(X;c,d,k),  if 

1 -'p 

♦ 

- — -)F.(h(x))-h'(x)f^(h(x))(^  - ^)F,(x)=0, 
1 j - - 'i  •'•'j  i 

for  i • j , i , j = 1 , . . . ,k. 


pjpP— (I.),:  “ F(x-':(0)>  - 'SO  " and  h(x)  = x+d,  d > 0 

where  r(\)  is  a Schur-convex  in  v , then 


(3.2.31)  F (x) 

•'S 

In  this  case,  (3.2.30) 


'SO  xSx) 

= -f(x-SO)  ---  - -f^(x) 

■ 'i  1 '■  i 

is  satisfied  and  i.(sc,d,k)  satisfies  (3.2.29). 
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Case  (1i):  If  F (x)  = ^ ^ ® 

c ji  1 where  n(A)  is  a Schur-convex  in  a,  then 


(3.2.32) 


3F^(x) 


'•'i 


^ n(A)]  = - f.(x) 

( ',  ] " ' i " n I A i A 


dn{A ) 


J A j 


In  this  case,  (3.2.30)  is  satisfied  and  v(A;c,d,k)  satisfies 
(3.2.29). 


Actually,  in  the  above  two  cases,  v(>;c,d,k)  does  not  involve 

A,  i.e.,  '^(x;c,d,k)  is  independent  of 

Let  S denote  the  size  of  the  subset  selected  by  the  procedure 

Rpi.  Let  .Q'  = - (a  1 , . . . , \j^) : ..  C E,  i = l,...,k  and  there  exist 

t,  such  that  A,  • a.  ^ ...  • , :. 

' 'Mni'"2ni  m '""k 

Let  E (S'R,  ) be  the  expected  size  of  the  selected  subset  using 
(jj  n 


'h’ 


Let 


-[1]  : -[?] 


-[k]’ 


for 


'.I  is  easy 


to  see  that 


(3.2.33) 
where 

(3.2.34) 


E (s:rj  = V p 

i = l 1 


p.  = i ::  F (h(x))dF,  (x), 

' i r-1  ^[r]  -[i] 

r*i 


and  F,  is  the  c.d.f.  of  T,.^  which  is  associated  with  Ar--), 

Aj-^-j  vU  “L'J 

i-l,...,k.  Using  the  same  arguments  as  in  Gupta  and  Panchapkesan  [37] 
and  Theorem  3.?.?,  we  have  the  following  theorem. 


Theorem  3.2.4.  Let  ■F.(x),  \ C F-  and  >F  (h(x)),  \ g E'  be  as  in 
the  hypothesis  of  Theorem  3.2.3.  For  . i.  F (S'R^)  is  Schur-convex 


i 


1 


s 

! 


'hL 
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in  X,  where  Xj-^^  =...=  x^^^ 


-[m+1]  ^ ^[k]  consequently 


sup  E(S|R^)  takes  place  at  V where  x^j  =...=  X|-j^^  = V provided  that, 

for  x-j  < X2  and  x-j,  Xp  e E,  the  following  holds; 
m 

(3.2.35)  (h(x))f.  (x)-(^  - v^-)F.  (x)f  (h(x))- 

dX.,.  dX^.  A.,  X^  dX.,.  oX.,.  X^  X2 

h'(x)  > 0, 


for  i < j,  i,j  = l,...,p  where  X-  = (x.^ ’^ip)’  ‘'  = '■’2. 


f,  (x) 
-1 


Remark  3.2.1 . (i)  If  ^ is  nondecreasing  in  x for  x.j  < X2,  then 


(3.2.35)  is  satisfied,  when  (1)  E;^(x)  = F(x-f;(x)),  -■»  < r,(x)  < «= 

and  h(x)  = x+d,  d > 0 where  f.(x)  is  Schur-concave  in  X,  or  when 
(2)  F,(x)  = F(-77-r)>  X > 0,  n(x)  > 0,  and  h(x)  = cx,  c > 1 where 
n(x)  is  Schur-concave  in  x. 

(ii)  If  (3.2.35)  is  satisfied,  then 


(3.2.36) 


sup  E (SjR, ) - k sup  P[CSjR.  ] 


where  = {/£  = (^1  ’ • • • ’^:k ) ■ t]  ^k’  -i  ^ 

(iii)  If  (1)  F.(x)  = F(x-f,(x)),  < f,(x)  < «>  and  h(x)  = x+d 

A ~ “ 

where  r,(x)  is  Schur-concave  in  a or  (2)  F^^(x)  = F(;:j-^y)%  x 0> 

n(x)  > 0 and  h(x)  = cx,  c 2.  ^ where  n(*)  is  Schur-concave  in  x, 

f;,  (x) 

-1 

also  if  is  nondecreasing  in  x for  x.j  .\2,  then 


(3.2.37) 


sup  E (SIR,  ) = k P*. 
, iij  ' n 


3.3.  A sufficient  condition  for  the  monotonicity  in  terms  of 
weak  majorization 

First  we  give  the  definition  of  weak  majorization. 

Definition  3. 3. i . A vector  a = is  said  to  weakly 

majorize  a vector  b = (b| , . . . ,bp) , if  a^  ^. . ap,  b^  ;>.  . .>  bp 
r r 

and  I a.  > } b • , r = l,...,p.  If  a weakly  majorize  b,  then 

i=l  ^ “ i=l  ^ 

we  write  a >>  b. 

m 

We  state  the  following  theorem  which  is  the  characterization  of 
weak  majorization  (see  Nevius,  Proschan  and  Sethuraman  [56]), 

Theorem  3.3.1.  Z >.■>  Z'  if  and  only  if  f(Z)  ^ f(Z')  for  all 
m ' 

nondecreasing  Schur-convex  functions,  where  f(Z)  is  defined  for 

Z,  1...1  z^. 

Let  X = (X,  ,...,X.  ) be  a random  vector  of  independent 
- n 'p 

components  where  x c E and  the  distribution  function  F,  (x)  of 

'i 

random  variable  X.^  is  stochastically  increasing  in  i = l,.... 
Let 

(3.3.1)  T = g(X^  ,...,X,^  ) be  a function  of  X,  ,...,X,^  . 

''1  'p  '1  'p 

Let 

(3.3.2)  F (x)  - P[T  - x]  - PTQ'X,  ) < xj. 

In  the  following  theorem,  we  obtain  a sufficient  condition  for 
E., (T,>)  to  be  nondecreas i ng  and  Schur-convex  in  \ c E for  some 
function  ,(T,')  as  defined  in  Thnoreis  :'.3.2. 
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’’neorem  3.3.?.  If  iP  (x) 


A ) e Ej  is  a 'airily 


continuous  distributions  on  the  real  line  such  that  F,(x)  is  a 
differentiable  function  in  x and  each  F,^  (x)  is  stochasti cal  ly 
increasing  in  and  is  a bounded  real-valued  and  differen- 

tiable function  in  x and  a^,  i = l,..,,p,  then  E(i,(T,A)  = Ev(g(X,),-) 
is  nondecreasing  and  Schur-convex  in  t E provided 

(3.3.3)  (i)  ’+;(g(X,),A)  is  nondecreasing  in  and  X,  and 

(ii)  c-tF,(x)(7-^ — — )^Hx,a)  - -77  '*(x,a)(^  - t-)^)F..  (x)>0, 

.X  A ,A.  aA.  - .X  - a...  a.,,  , 

for  i < j , i ,j  = 1 , . . . ,p. 

Proof.  3y  Theorem  3. 2. 2. and  (3. 3. 3)-(i  i ) , E;;;(T,x)  is  Schur-convex 


Now  we  want  to  show  that  E.'(T,x)  is  nondecreasing  in  \.  Let 


. _ f • 

V ‘ 1 1 


P'  ’ I 


Assume  X > x'.  We  define  • 


; ' 'M A^),  >1  i...>  ,A. 

• ■ •'  in  the  sense  that  x^  2^  i = 


l,...,p.  Since  F.  (x)  associated  with  the  random  variable  is 
' i ' i 

stochastically  increasing  in  i.e.  if  \,.  7 i = l,...,p, 


(3.3.4) 


i St  T 


. i - 1 ,P. 


Since  X,  ,...,X.^  are  independent  random  variables,  then 
' 1 P 

(3.3.5)  X,  • X,  . 


By  (3.3.3)-(i),  since  ';(g(Xj,:.;  is  nondecreasing  in  X_,  then 


(3.3.6' 


.(q(X.),0  c fg(x, ,),•). 


Since  ij;(g(X  . ) ,a  ' ) is  nondecreasing  in  x ' and  x ^ X ' , 
(3.3.7)  ii»(g(x^.),x)  > i|/(g(x^,),x'). 


From  (3.3.6)  and  (3.3.7),  we  have 

(3.3.8)  <j.(g(Xj,x)  > ^(g(x.  ,),x'). 

■-  ■ St  ■ 

Hence 

(3.3.9)  E,j,(T,x)  > E4,(T.x‘). 

This  completes  the  proof. 

By  applying  Theorem  3.3.2,  we  can  obtain  the  following  results. 

Corollary  3.3.1.  If  F (x)  is  Schur-concave  in  x,  <^(x,x)  satisfies 

A " * 

(3.2.26),  t(;(x,x)  is  nondecreasing  in  x,  v(g(X^),x)  is  nondecreasing  in 
and  X and  F^  (x)  is  stochastically  increasing  in  X^,  i=l,...,p, 

then  Ei|/(T,x)  is  nondecreasing  and  Schur-convex  in  X. 

Corollary  3.3.2.  If  F (x)  is  Schur-concave  in  x,  ij;(x)  is  nondecreasing 

A ”* 

in  X,  4<(g(X^))  is  nondecreasing  in  X^  and  F^^  (x)  is  stochastically  increasing 

in  X.,  i=l,...,p,  then  Ei|i(T)  is  nondecreasing  and  Schur-convex  in  x. 

k-1 

By  letting  j,(x,x)  = F (h(x))  and  using  Theorem  3.3.1  and  Theorem 

” A 

3.3.2,  we  have  the  following  theorem. 

Theorem  3.3.3.  Let  {F,(x),  a C E}  be  a family  of  continuous 

A "" 

distributions  on  the  real  line  such  that  F (x)  and  F (h(x))  are 

A A 

differentiable  functions  in  x and  x^.  and  F^  (x)  is  stochastically  increasing 
in  x^,  i = l,...,p.  Let  ij)(>;c,d,k)  be  defined  in  (3.2.10). 

Then  4)(x;c,d,k)  is  nondecreasing  in  x C E in  the  sense  of  weak 
majorization , i.e.,  it  is  nondecreasing  and  Schur-convex  in  x, 
provided  that 
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(3.3.10)  (i)  F^(h(g(Xj^))  is  nondecreasing  in  and  x and 

- 3X7^ 

**  i 3 * ^ J 

F^(x)  ^ 0,  for 

i < j.  i.j  = 1 P,  where  fj^(x)  = ^ 

3.4.  Selection  of  the  population  associated  with 

If  the  best  population  is  defined  to  be  the  one  associated 
with  Xp,-i,  where  Xr,-i  < x^ , i = 1 ,. . . ,k.  We  now  define  a class 

of  procedures  for  the  selection  of  the  population  associated 
with  x^^ip 

Let  H = c € [1,“),  d € [0,«)  be  a function  defined  on 

the  real  line  satisfying  the  following  conditions.  For  every  x 


(3.4.1)  (i)  H^^^(x)  < X 

(ii)  H^^q(x)  = x 

(iii)  j(x)  is  continuous  in  c and  d 
(iv)  j(x)  + as  d ->  “>  and/or  xH^  ^j(x)  + 0 as  c 

Of  particular  interest  are  the  functions  x-d  and  - d. 

A class  of  procedures  for  selecting  a subset  containing  the 
best  is  defined  as  follows. 


R^;  Select  population  if  and  only  if 


(3.4.2) 


H(T.. ) £ min  T . 
^ l<r<k  ^ 


The  probability  of  a correct  selection  is  given  by 
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« k 

(3.4.3)  P[CS|RJ  = / n r (H(x))dF,  (x). 

^ - r=2  ^[r]  ^[1] 

where  P^(x)  = l-F^(x). 

Because  of  the  assumption  (3.2.8), 


(3.4.4)  P[CS|R.]  > /”  (H(x))dF  (x). 

H ^[1]  ^[1] 

Let  n = {o)  = (Xi,...,X(^):  a e E,  i=l,...,k  and  there  exists  some 

” ~i 

i such  that  A^  < A-  V^}. 

Hence 

(3.4.5)  inf  P[CS1R.]  = inf  <p(A;c,d,k) 

n A ■ 

where 

(3.4.6)  cp(A;c,d,k)  = / (H(x))dFj^(x)  and  A € E. 


Using  the  same  method  of  proof  as  in  the  case  of  R^,  we  have  the 
following  results. 


Theorem  3.4.1.  For  the  procedure  R^^,  cp(A;c,d,k)  is  Schur-convex 
in  A € E,  provided  that 

(3.4.7)  (^  - ^)F^(x)-H'(x)f^(H(x))-f^(x)(^  - g|^)F^(H(x) )>0, 

for  i < j , i , j=l ,. . . ,p,  where  H' (x)  = ^ H(x). 

Let  n'  = {w  = (Ap...,A|^):  A.  6 E,  i=l,...,k  and  there  exist 

ii, such  that  A > A > ...  > A }. 

' ^1  m *^2  m m “^k 

Theorem  3.4.2.  sup  E(S|R^)  takes  place  at  V where  A^  =...=  Aj^  = V 

provided  that,  for  A-|  < A2»  a-j,A2  G E, 

m 
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(3.4.8)  H'(x)f^^(H(x))(^  - 


F,  (H(x))  > 0, 
^1 


for  i < j,  i,  j=l,...,p  where  x-  = ■ .^,-p) . i = 1,2. 

3.5.  Selection  procedures  for  multivariate  normal  distributions 
in  terms  of  majorization. 

Let  TTp...,ir|^  be  k populations.  Let  it^  be  associated  with 
-i  “ (X-n  •• • • *Xjp) > i = l,...,k  where  is  random  vector  with  a 
p-variate  normal  distribution  with  unknown  mean  vectors  = 

(u.j^ . ,ji.jp)  and  positive  definite  covariance  matrix  l.  We  assume 
that  V ^ ^ 1 ""  l,...,k. 

It  is  assumed  that  among  the  k given  populations,  there  always 
exists  u^i^j  such  that  for  all  i = l,...,k.  This  popula- 

tion is  called  the  best  population. 


Let  denote  n independent  observation  vectors, 

each  with  p components,  from  population  it.. 

Let 

(3.5.1)  X^  - (X^-jj^,  ^-j  2p  ’ " * " ’ ^ipA^  ’ ^ ” l,-**,k,  A - l,...,n. 
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[ 

1 


L 


(3.5.2) 

1 n 

*“  “ / ^ l>**»»k*  j — l«***fP» 

(3.5.3) 

''ij 

- ^ > i ~ l»-*.*k>  j ~ l(...iP> 

«,=!  ^ 

(3.5.4) 

s?. 

U 

" n U^ijm  ■ " ■'•••••k.  j = 1.....P. 

m=  1 

(3.5.5) 

(3.5.6) 

~ +•-.+  ^ ~ j — 1y«»*9p9 

(3.5.7) 

S[1] 

w(i) 

j 

~ ^(0^  ^ ~ j *“  1»»***P* 

Let 

us  denote  by  the  population  (unknown)  associated 

with  mean  vector  and  let  observation  associated 

with  i = l,...,k,  j = Our  goal  is  to  select  a subset 

(may  be  empty)  of  the  k populations  so  as  to  include  the  popula- 
tion associated  with 

(A)  Assume  i is  known. 

We  propose  a selection  rule  as  follows: 

R-| : Select  population  if  and  only  if 

(3.5.8)  Y..  > max  Y . - d,  j = l,...,p-l  and 

~ l<v<k 

min  Y _ + d > Y.  > max  Y -d. 

l<v<k  ''P  - IP  - i<^<k 

Let 

(3.5.9)  Q = D E*  D' 


((k-l)p)x(kp) 


A 0 B 


0 0 0 A B 


l* 

(kp)x(kp) 


T 0 


0 T 


. T = ^ CEC 
pxp 


■ "V  ■■  (■'"  1 

) * w pxp  \o  ‘-7 


and  C = I 1 1 

V : ‘m 

\i  1 1 


Let  fi  = {u  = (ui €E-|,  i = l,...,k  and  there  exists 


i such  that  u.  > V j),  where 


(3.5>10)  E^i  - {A  ~ (Aij...jAp)j  X"!  ^p*  ^ 

The  following  theorem  gives  the  infimum  of  the  probability 
of  a correct  selection. 


Theorem  3.5.1 . 

(3.5.11)  inf  P[CS|R,]  = P[Z.<d,|Z.(<  d,  for  i e {1 ,2,. . . ,(k-l)p}- 
fi  ' ' J 

{p,2p,. . . ,(k-l  )p}  and  j e {p,2p,... ,(k-l)p}]. 


where  i • • • .Z(|^_1  jp) ' - N(  0 , Q ). 

(k-l)pxl  (k-1 )px(k-l )p 
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Proof. 


(3.5.12)  P[CS|R^]  = P{Y(,^)j  I Y(|()p-Y(^jp|  < d,j=1 ,. . . ,p-1 , 


Y.j  - N(vj,  i = l,...,k,  j=l,...,p,  where 

~ (1»-**»1>0,...,0),  j = l,...,p. 
j times 

Since  i = l,...,k-l, 

(3.5.13)  inf  P[CS|R,]  = P{Y*.  > Y*.-d. [ Y*  -Y*  |<d,  j=l,....p-l. 


where 

(3.5.14) 


(Y*i....,Y*  )•  - N(0,T),  V = 1 k. 


V“ 1 9 ■ • • y k“ 1 } 9 


!:=  •■■.’'Ip)'.  ” ■ '.■■■.!<■ 

Hence  from  (3.5.14), 

(3.5.15)  Y = (Yf,...,Yj^)'  -•  N(  0 . E*). 

kpxl 


Let  = Y*.  - Y*^. , V - l,...,k-l,  i - 1,...,p. 


^^vl ' ''  l>..,k-l 


Since 


Z = (Zt, . . . ,Zf  1 ) ' = D Y and  from  (3.5.15),  then  Z - N(0,Q). 
(k-l)pxl 

This  completes  the  proof. 


Theorem  3.5.2.  If  Q = De*D'=  (q^j)  is  positive  definite  with 

2 2 

‘^11  "•••■  *^(k-l  )p,(k-l )p  ""  " '^ij  = O'  » i f j,  a and  p are 
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known,  then 

(3.5.16)  sup  E(S|R, ) = KP*  provided  that  inf  P[CS|R,]  = P*. 

n n 

Proof. 

(3.5.17)  E(S|R,)  = P(V,^  > -d.  Y^p-d  < < Y^^  Y d. 

V^i  $ ^ ” l9.«.tlCs  » • • • »p“l  } 

' J,  i"- 

Where  - Y.^)  - (V^-vj).  vfi.  v.i=l k.  t=l p. 


(3.5.18)  Z.  = {Z.^^}  - N(0,Q). 
(k-l)pxl 


B.j  “ ^--j'  — »*»*»P“1»v^i,v  “ l,...,k}. 


li  ^ivJi  " K'K^  v=^i,v,  i=l,..,k,a=l,...,p. 

lx(k-l )p 

Hence 


(3.5.19)  E(S|R^)  = P[B.+T^.]. 


Since  the  joint  density  of  Z which  is  defined  in  (3.5.18)  is 

Schur-concave  (see  Marshall  and  Olkin  [50  ])  and  since  ^ £ B.|  and 

X < y implies  x € B. , then  by  Theorem  2.1  of  Marshall  and  Olkin  [50], 
' m ~ ' 

P[B.+I^. ] is  Schur-concave  function  of  T..,  i = l,...,k. 
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s' 


i 

\ 


5 


i' 


Since  T.  > (a^,...,a^),  i = where 

m 


A^.  = vj  + ...+Vp _i,  then  from  (3.5.19) 
k 

(3.5.20)  E(SlR^)  1 P[B^+(a.,...,a.)] 


where  a^j^-j  is  the  largest  value  among  {a^ 
sup  E(S|Ri)  is  obtained  when  a^  = ...=aj^. 
=...=  A|^.  Hence  a^  =...=  aj^  = 0. 


If  a^  =.. 


That  is, 

= aj^,  then 


Thus 

k 

sup  E(S|Rj  = I P[BJ  = k P*  provided  that  inf  P[CS|R,]  = P*. 
' i=l  ’ n 

This  completes  the  proof. 


2 0 

(B)  Let  z = , n)  where  o is  known.  Without  loss  of  generality, 

0 * 

we  can  assume  a = 1.  We  propose  a subset  selection  rule  as 
follows. 

R2:  Select  TT^.  if,  and  only  if, 

(3.5.21)  Y..  > max  Y . - for  j = l,...,p-l  and 

~ l<v<k  /n 


V 


I 
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Since  vj"'  i Vj''^  i = l,...,k-l,  then 

infPtcsiR^]  -p[z„,  id.iz^plid,  ] 

where  Z^^-,  v = l,...,k-l,  i=l,...,p  are  N(0,1)  random  variables 
with  covariances  as  defined  in  (3.5.23).  This  completes  the  proof. 


If  k = 2 and  p = 2,  we  can  show  that 


(3.5.26)  inf  P[CSlR2]  = / 


(1-/1)^^ 


)-$(- 


-d-(y) 


I\l/4 


d-(i)'/^x 


(1-/J)^  (1-/^)^' 


-)]di>(x) 


Theorem  3.5.4. 

(3.5.27)  PECSIR^]  > max{[4.(d)-«(-d)]P^'^'^\l-(p+l)(k-l)<!.(-d)). 


Proof.  Let  {Z  •}  be  defined  as  in  Theorem  3.5.3,  then 

Vl 

(3.5.28)  P[CS|R2]  i P[Z^.  5 d,!Z^pi  < d,  i=l p-1 ,v=l , . . . ,k-l ] . 

Note  that 

m mm 

(3.5.29)  P[nA.]>l-y  zl  P[B.  J 

i=l  ^ i=l  j=l 

where  A,,..., A denotes  a seoo'r'ice  of  events,  and  A^I  the  event 

n 

complementary  to  A.,  such  that  = U , i=l , . . . ,m. 

1 1 j=l 

Hence  from  (3.5.28)  and  (3.5.29), 

k-1 

(3.5.30)  P[CS|R2J  i 1-  I P[Z^i>d]-  I P[]Z  J>d] 

v,i  v=l  ^ 

= l-(p-l)(k-l)^(-d)-(k-l)2$(-d) 


- l-(p+l)(k-l).i.(-d). 
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( 

i 

I 

[ 


By  [ 68  ] and  from  (3.5.28) 

(3.5.31)  P[CS|R2]  lP[|Z^.j|<d.  v=l , . . . ,k-l  ,i=l . . . . ,p] 

> ^ P[|Z^.l<d]  = [$(d)-$(-d)]P(''-''^ 

v,i 

From  (3.5.30)  and  (3.5.31),  we  prove  the  theorem. 

For  given  P*(^  < P*  < 1),  a conservative  value  of  d (d  > 0) 
can  be  obtained  by  letting 

max{[$(d)-4.(-d)]P^'^‘^  ^ l-(p+l)(k-l)$(-d)l  = P*. 

Let  d.|  > 0 be  the  value  such  that  [4>(d)-4>(-d)]P^*^'^  ^ = P*  and 
let  d2  > 0 be  the  value  such  that  1-(p+l  )(k-l  )4'(-d)  = P*,  then 
the  minimum  d(>  0)  satisfying  the  basic  requirement  is  given  by 
d = min{di ,d2) . 

Theorem  3.5.5. 

(3.5.32)  E(SlR2)  < k[$(d)-<f.(-d)] 


Proof. 


Y.-,-  > Y , - 


d/ZT 


01  - VI  ^ 

I 

(3.5.33)  E(S|R2)=  I P 

* y _ < Y . < Y + 

vp  ^ - jp  - vp  ^ • 


for  i=l , . . . ,p-l j ,v=l , . . . ,k  }, 


Y .-Y  • •-(VV,-V'?) 

Let  Z..  i - l,...,p,.:tj,.=l k. 


'Vl 


/?7//n 


Then 
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k V-^-v'^ 

(3.5.34)  E(S|Rp)  = y p-z  . • d + ,IZ  J<d, 

^ j'^1  “ mi/^  "P'- 

1 = 1 ,. . . ,p-l  ,vij,v=l  ,. . . ,k} 


1 X P^^ui  1 ^ •'  7 ■^0'"  ''+j  .v=l ,. . . ,k}  = A.  (say), 

.1  = 1 v'ly  v'n 


1 = 1 p-1 . 


A1  so. 


E(SIR2)  < _y^P[iZ,^p:<d,v4j,v=l,2,...,k]  - A (say). 

vJ-vV 

Since  P[Z.<d  + v4j , v=l , . . . ,k] 

1 

1 z P^Z.  • < d -1 ], 

' v=l  - z?i7/iT 

then  vfj 

k k V'^-V'( 

(3.5.35)  A.  < I / (d+ — ^ — ^)  (since  Z . ' N(0,1)) 

^ ' j = l v)=l  /FT//n 


" ' j = l V=1  m//n 


= • Q (say) 


1 k 

where  order  values  of  {V  ■ , . . . ,V^  . 

By  the  same  argument  in  Gupta  and  Huang  L34],  we  can  show 
that  the  sup  of  Q is  obtained  when  Ap.|j^  =...= 

Thus  from  (3.5.35) 

k k 

(3.5.36)  A.  - pL  y ) .!=(d)  = k4.(d),  i = l p-1. 

j=l  V=1 


Also 


(3.5.37)  Ai  PCIZ^pl  < d] 


= k[t>(d)-<t(-d)]. 

Hence  from  (3.5.36)  and  (3.5.37), 

(3.5.38)  E(S|R2)  1 min{k<H'd),  k[4>(d)-4>(-d)]} 

= k{$(d)-<t>(-d)}. 

This  completes  the  proof. 

2 0 2 

(C)  Assume  i = (°  • , ^,)  where  is  unknown. 

0 • 

We  propose  a subset  selection  rule  as  follows. 


R^:  Select  population  if  and  only  if 

^ d R 

(3.5.39)  Y..  ^ max  Y . j = l,..,,p-l  and 


^2  d R.  /2  d R. 

max  Y ^ ^ < Y.  < min  Y „ + 

V /fTT  " 'p  - V /frr 


Theorem  3.5.6. 


(3.5.40)  P[CS|R.,]  > l-(k-l)  [ p[t.  < -d]-(k-l)P[|tJ  > d] 

-5  i=l  ’ P 

where  t..  is  a r.v.  having  a t-distribution  with  i(n-l)d.f.,  i=l. 


Proof. 


(3.5.41)  PCCSIRj]  . - 


JZ  d 


’(k)P  - ’Wp  * ^''' 

v=1,...,k-l] 

= P[t  . < d - -J ^ < d,j=l P-I.v=l,...,k-1] 

■^"OOJ  “ 


where 


(3.5,42)  t„j= 


W^-<''‘'’’-''<''’' 


1 / 

• >^n“l  y j“l,,..>p*V"l»...*k-l. 


5’"PP  <''(v)j-''(k)j>  - N(V<">-V<^).?i,2, 

P^{|/\i  2 

and  — is  distributed  as  a r.v.  of  x with  j(n-l)  d.f, , then  t . 

o ^ “ 

is  a Student's  t-distribution  with  j(n-l)  d.f. 

Since  = l,...,k-l,  j = l,...,p,  then 

(3.5.43)  inf  P[CS1R3]  = ?[t^.  i<l,!t^pl  id.  j=l , . . . ,p-l  ,v=l k-1] 

1 1-  I P[t  i 1 d]-  [ P[|t|<  d]  (by  (3.5.29)) 

v,j  \>=] 

P-1 

= l-(k-l)  I P[t.  < d]-(k-l)P[|t_|  < d] 
j=l  P 

where  t-  is  a r.v.  of  a t-distribution  with  j(n-l)  d.f.,  j=l,...,p. 

This  completes  the  proof. 


Theorem  3.5.7. 


(3.5.44)  E(S|R3)  <min{kG^(d),  k[G  (d)-G  (-d)]} 


where  Gj(x)  is  the  c.d.f.  of  a r.v.  of  a t-distribution  with  j(n-l) 
d.f.,  j = l,...,p. 


no 


Proof. 


dv^  R.. 
Y..  > Y , ^ 


1J  - vj 


/n-1 


(3.4.45)  E(S|R3)=  P 


R. 


^ < Y._  < Y 


vP  ^ - ■'P  - ^P  /fTT  ’ 

for  j = l,...,p-l,  vfi,  \)  = 1 k}. 

Y ,-Y.  .-(v'^-V^. ) 

Let  t . = — ^ • /n-1 , then  t . is  a Student's  t- 

,^R.. 

^ J 

distribution  with  j(n-l)  d.f. 

vV-vl 


(3.5.46)  E(SiR.)=  I P 
i = l 


Id. 


Hence 


v!-v^ 


(3.5.47)  E(S|R,)  < \ P{t  . < d + -J— ^ • /h-'l , vfi,  v=1,...,k} 

' 3 - i^i  vj  - ^ D 


1.1 

Aj  (say),  j = 1.....P-1. 


and 


(3.5.48)  E(S|R^)<  I P{|t  1 < d,  vfi,  v=l,...,k} 

i i=T  vp 

= A (say). 

As  in  Theorem  3.5.5,  we  can  show  that  the  sup  of  A.  is  obtained 

V 

when  VI  =...=  v![,  for  j = l,...,p-l.  From  (3.5.47)  and  (3.5.48), 

J J 

we  get  A.  < kG.(d),  j=1,...,p-1  and  a ktG„(d)-G_(-d)].  Thus 

J J P P 

E(S1R3)  < minfk[Gp(d)-6p(-d)],  kGj(d),  j=l p-1) 

= min{k[Gp(d)-Gp(-d)],  kG,(d)} 


in 


3.6.  Selection  procedures  for  multivariate  normal  distributions  1n 
terms  of  weak  majorization 

It  is  assumed  that  among  the  k populations,  there  always 
exists  such  that  » Uf,].  i = l,...,k.  This  population 

is  called  the  best  population.  We  are  using  the  same  notation  as 

in  Section  3.5.  Let  = {w  6 (wi ,. . . .Pj,,):  p^-  6 E-j,  i=l,...,k  and 

there  exists  some  j such  that  p^  » p<  Vi}  where  E,  is  defined 

m " 

in  (3.5.10). 


(A)  Assume  r is  known. 


We  propose  a selection  rule  as  follows. 

Rg:  Select  population  if  and  only  if 

(3.6.1.)  Y..  > max  Y .-d,  j = l,...,p. 

’J  l<v<k 

Using  the  same  argument  as  in  Theorem  3.5.1,  we  obtain  the  following 
result. 

Theorem  3.6.1. 

(3.6.2)  inf  P[CS|R.]  = P[Z.  < d,  i = 1 ,. . . , (k-1 )p] 

3 1 


(Zi ,. . . ,Z(|^_^  jp) ' - N(  0 ,Q),  where  Q is  defined  in  (3.5.9). 

(k-1 )pxl 

2 Q 

(B)  Assume  E * (°*,  ^,)  where  a is  known  and  we  assume  o = 1. 

0 'c 

We  propose  a selection  rule  as  follows. 

Rgi  Select  population  ir^  if  and  only  if. 


(3.6.3.)  Y..  > max  Y ^ , j«1,...,p. 

l<v<k  ''J  ^ 


It  is  similar  as  in  Theorem  3.5.3,  we  can  show  that 


Theorem  3.6.2. 

(3.6.4)  inf  P[CS|RJ  = P[Z  . < d.  v=l....,k-l,j=l,...,p] 

where  is  defined  in  Theorem  3.5.3. 

If  k = 2,  p * 2,  we  can  show  that 

d_(l)l/4j^ 

(3.6.5)  inf  P[CS|R,]  = ]il»(x). 

" (1-/  J)* 

In  a same  manner  as  in  Theorem  3.5.4  and  Theorem  3.5.5, 
we  have  the  following  results. 

Theorem  3.6.3. 

(3.6.6)  P[CS|Rg]  > max{[4.(d)-$(-d)3P^'^‘''\  l-p(k-l)[l-.>(d)]}. 


Theorem  3.6.4. 


(3.6.7) 


E(S|Rg)  < k$(d), 


2 0 2 

(C)  Assume  E = (°»,  %)  where  is  unknown. 

0 *0 

We  propose  a subset  selection  rule  as  follows. 
R,:  Select  population  it.  if  and  only  if 

(3.6.8)  Y^j  > max  Y^j  - ^ '^ij  , j = 1 


In  a manner  similar  to  that  in  the  proof  of  Theorem  3.5.6, 


we  have 
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